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Preface 


The  formal  goal  of  this  study  was  to  add  a  model  of  a  secondary  coolant  loop 
to  an  existing  model  of  a  space  reactor.  The  need  for  this  model  is  to  simulate  the 
operation  of  the  SP-100  space  reactor  and  to  confirm  the  preformances  reported 
by  General  Electric. 

This  study  is  an  extension  of  the  work  presented  by  Jong  T.  Seo.  Elements 
of  his  thesis  work  have  been  reproduced  in  this  study.  Most  notably,  I  presented 
descriptions  of  the  pump,  heat  exchanger,  and  thermoelectric  generator  models  in 
much  the  same  fashion  as  Seo. 


The  design  of  the  SP-100  has  changed  since  Seo’s  work,  it  was  undergoing 
changes  during  my  work,  and  will  likely  be  changed  again  in  the  future.  The  util¬ 
ity  of  this  study  is,  therefore,  appreciated  by  those  who  will  make  continued  stu¬ 
dies  of  the  SP-100  as  its  design  evolves. 


My  appreciation  and  thanks  go  to  my  faculty  advisor,  Lt  Col  R.F.  Tuttle, 
for  his  guidance.  Also,  thanks  go  to  Maj  JA.  Lupo  for  his  computer  expertise,  to 
Lt  Col  W.P.  Baker  for  help  in  matters  of  mathematics,  and  to  LCDR  K.A. 
Mathews  for  his  explanations  of  reactor  systems.  I  extend  my  compliments  to 
Jong  Seo  for  his  work  on  SNPSAM  and  thank  him  for  his  discussion.  A  word  of 


thanks  is  owed  to  Keith  Hoffman,  at  Foreign  Technology  Division,  for  providing 
me  with  an  unclassified  Soviet  perspective  of  space  reactor  technology.  Special 
thanks  goes  to  the  staff  of  AWYS  at  the  Air  Force  Weapons  Lab  for  sponsoring 


my  work  and  providing  assistance  and  to  the  SP-100  office  of  General  Electric  for  _ 

For 

their  interest  and  help.  Finally,  my  thanks  to  professor  C.A.  Sparrow  for  his  help  T 
and  selection  of  this  project.  d 


Vaughn  H.  Standley 
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Abstract 

The  purpose  of  this  study  was  to  add  a  secondary  coolant  loop  model  to  the 
Space  Nuclear  Power  Systems  Analysis  Model  (SNPSAM).  The  heat  rejection  sys¬ 
tems,  including  the  TEM  pump,  energy  conversion  assembly,  and  radiator,  are 
emphasized  while  the  reactor  model  is  de-emphasized.  The  specific  configuration 
chosen  for  the  secondary  is  based  on  literature  obtained  from  General  Electric. 

Modifications  were  made  to  the  TEM  pump  and  heat  exchanger  models.  A 
subroutine  that  simulated  the  secondary  was  written  and  SNPSAM  was  modified 
to  incorporate  this  subroutine.  Modifications  included  subroutines  for  computing 
head  loss,  heat  transfer  coefficients,  and  mass  transfer.  In  the  revised  SNPSAM, 
the  steady-state  mode  generated  rejected  power  outputs  that  defied  an  energy  bal¬ 
ance  in  the  system.  The  cause  for  the  imbalance  is  in  the  iterative  method  used 
by  the  thermocouple  model.  The  iterative  method  used  in  the  thermocouple 
model  made  SNPSAM  sensitive  to  initial  conditions.  In  the  transient  mode  opera¬ 
tion,  the  model  generated  solutions  that  were  accurate  but  displayed  oscillations 
in  the  first  few  seconds  of  simulation.  The  steady-state  performance  predicted  by 
SNPSAM  was  compared  with  the  steady-state  performance  reported  by  General 
Electric.  Since  the  steady-state  mode  would  not  function  satisfactorily,  a  slowly 
varying  transient  case  was  used  in  its  place.  Emphasis  is  given  to  procedures  that 
facilitate  tailoring  the  model  to  whatever  design  is  current. 
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CUSTOMIZING  SNPSAM: 
INTRODUCING  A  SECONDARY  COOLANT  LOOP 

I.  Introduction 


1.1  Background 

Military  and  civilian  efforts  in  space  will,  in  the  near  future,  place  demands 
on  auxiliary  power  systems  that  conventional  designs  cannot  meet.  Greater  power 
and  longer  life  are  two  criteria  that  challenge  future  systems.  A  sensible  alterna¬ 
tive  is  nuclear  power.  The  necessary  ingredient  for  higher  power  is  the  develop¬ 
ment  of  a  system  that  will  provide  an  energy  density  in  the  range  of  1,000  to 
3,000  k\Vc /kg  (12:13).  Nuclear  reactors  lend  themselves  to  this  range  of  power 
densities.  Under  a  Department  of  Defense  contract,  General  Electric  proposes  to 
design  the  SP-100  portable  space  nuclear  reactor.  Figure  1.1  illustrates  the 
configuration  of  the  SP-100.  The  SP-100  is  a  passive  reactor  design.  With  few 
moving  parts,  the  SP-100  promises  long  life  and  predictable  operation. 

The  SP-100  is  a  fast  neutron  spectrum,  highly  enriched  (93%  U-235),  liquid 
metal  reactor  (1:234).  The  core  is  composed  of  uranium  nitride  fuel  pellets  and 
niobium- 1  zirconium  (1  %  zirconium)  alloy  clad.  The  reactor  is  controlled  by 
drums  that  are  rotated  to  expose  more  reflector  material  {Be  0)  or  more  absorbing 
material  ( BAC ).  The  SP-100  reactor  will  be  actively  cooled  by  a  liquid  metal  cir¬ 
culated  by  thermoelectric  electromagnetic  (TEM)  pumps.  The  liquid  metal  will 
likely  be  lithium  since  it  has  high  electrical  and  thermal  conductivity,  has  a  low 
melting  point  and  high  specific  heat  capacity.  Heat  is  rejected  by  radiator  panels 
that  are  constructed  of  11,220  heat  pipes  welded  to  a  central  channel  through 
which  the  coolant  passes.  The  radiator  panels  are  deployable  to  increase  the 
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effective  surface  area  of  the  radiator.  Thermoelectric  generate  electricity  in  the 
SP-100.  The  thermoelectric  converters,  or  thermocouples,  are  fixed  between  the 
reactor  coolant  fluid  {the  hot  leg)  and  the  radiator  panels  (the  cold  leg). 

The  mission  of  the  SP-100  is  to  demonstrate  that  space  nuclear  reactors  are 
feasible  and  practical.  The  system  will  be  flown  by  the  shuttle  in  the  1990’s  and 
has  a  planned  lifetime  of  seven  years.  The  SP-100,  if  successful,  will  be  the  first 
U.S.  reactor  in  space  since  the  SNAP-10A  which  was  launched  in  1965*. 

The  SP-100  is  still  in  the  design  phase.  Dimensions,  configurations,  and 
materials  have  not  been  formally  accepted.  General  Electric  is  studying  various 
trade-ofis  between  system  mass,  power  output,  efficiency,  compactness,  scalability, 
longer  life,  and  cost. 

The  Air  Force  is  interested  in  conducting  its  own  study  of  the  SP-100.  As  a 
part  of  this  "in-house"  effort,  the  Air  Force  has  procured  a  computer  model  of  the 
SP-100.  The  model  is  called  Space  Nuclear  Power  Systems  Analysis  Model ,  or 
SNPSAM.  SNPSAM  was  written  by  Jong  T.  Seo  for  the  Institute  of  Space 
Nuclear  Power,  The  University  of  New  Mexico,  Albuquerque,  NM  under  contract 
with  the  Air  Force  Weapons  Laboratory. 

SNPSAM  was  written  to  simulate  the  steady-state  and  transient  operation  of 

the  SP-100  and  similar  space  reactor  systems.  SNPSAM  contains  4986  lines  of 

Fortran-77  code  and  reads  an  input  file,  entitled  SYS.INP,  that  is  167  lines  in 

length.  SNPSAM  writes  30  output  parameters  to  PLOT.TMP  for  each  time-  step 

in  a  simulation.  Samples  of  SYS.INP  and  PLOT.TMP  are  included  in  Appendices 

A  and  B,  respectively.  SNPSAM  was  written  to  model  a  specific  space  reactor 

configuration.  Changes  can  be  made,  however,  to  parameters  such  as  dimensions, 

materials,  and  environments (e.g.,  external  heat,  ambient  space  temperature,  etc.), 

*  The  SNAP-1QA  was  shut  down  43  days  after  start-up  because  of  a  faulty  voltage  regula¬ 
tor  in  the  Agena  boost  vehicle.  The  program  was,  nevertheless,  judged  a  success. 
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by  modifying  the  input  file. 

The  actual  Fortran  source  code  of  SNPSAM  was  written  by  the  DSNP, 
Dynamic  Simulation  of  Nuclear  Power,  computer  code  p re- compiler.  The  reactor 
model  comes  from  the  DSNP  library  (13:1-3).  The  other  models  used  in  SNPSAM 
were  written  by  Seo  and  then  submitted  to  DSNP  for  compiling.  DSNP  organizes 
the  models  according  to  a  schematic  entered  by  the  user  in  the  form  of  a  com¬ 
mand  file.  SNPSAM  was  written  for  the  VAX/785  computer  with  the  VMS  4.1 
operating  system  but  was  compiled  and  executed  successfully  on  the  Air  Force 
Institute  of  Technology’s  Classroom  Support  Computer  (CSC),  which  is  a  VAX 
11/785  computer  with  the  VMS  5.1  operating  system. 

SNPSAM  was  initially  written  to  emulate  the  SP-100  configuration  depicted 
in  Figure  1.2.  On  the  surface,  the  configuration  shown  in  Figure  1.1  is  not  radi¬ 
cally  different  to  the  configuration  shown  in  Figure  1.2.  The  heat  systems,  how¬ 
ever,  are  considerably  different. 

The  previously  proposed  design  for  the  SP-100  space  reactor  utilizes  heat 
pipes  to  remove  heat  from  the  primary  coolant  loop.  The  most  recent  proposal  by 
General  Electric  includes  a  secondary  lithium  coolant  loop  (7:29).  The  added 
secondary  is  typical  because  coolant  flows  in  close  proximity  with,  and  counter 
current  to,  the  flow  of  primary  coolant  while  in  the  heat  exchanger.  Furthermo.  e, 
the  primary  and  secondary  are  driven  by  common  pumps.  SNPSAM  does  not 
include  a  model  of  the  secondary  lithium  loop  proposed  by  General  Electric,  but 
models  the  previously  proposed  SP-100  design  where  heat  pipes  transport  waste 
heat  directly  from  the  primary  coolant  loop  to  outer-space. 

There  have  been  other  changes  to  the  design  of  the  SP-100  since  SNPSAM 
was  written.  The  reactor  has  been  scaled  to  generate  greater  thermal  power.  The 
size,  shape,  and  operation  of  the  radiator  panels  has  been  modified.  General  Elec¬ 
tric  proposes  to  make  the  panels  larger  and  more  deployable.  The  heat 
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exchangers  have  been  redesigned  and  there  are  four  times  as  many  pumps  in  the 
new  design. 

There  are  three  reasons  for  adding  a  secondary  coolant  loop  to  the  SP-100. 
The  first  arises  because  heat  pipe  technology  is  limited.  Many  long,  very  small  in 
diameter  heat  pipes  would  be  required.  The  efficiency  and  reliability  gained  by 
using  heat  pipes  is  outweighed  by  their  cost  in  terms  of  mass  and  dollars  (9).  A 
second  reason  for  adding  the  secondary  coolant  loop  is  that  it  allows  the  energy 
conversion  assembly  to  be  easily  separated  from  the  radiator  (i.e.,  a  more  modular 
design).  The  third  reason  is  that  in  a  counter-current  heat  exchanger  there  is 
nearly  a  constant  temperature  gradient  across  the  thermocouples  along  the  length 
of  the  exchanger  and  this  may  cause  the  thermocouples  to  operate  at  greater 
efficiency. 

1.2  Scope 

This  study  is  devoted  to  modeling  the  SP-100.  Since  an  exact  design  for  the 
SP-100  has  not  been  selected,  one  of  the  latest  proposed  designs  is  modeled.  The 
design  of  choice  is  detailed  in  (7:1-13).  Work  towards  this  end  was  carried  out 
with  the  following  objectives: 

1.  Add  a  secondary  coolant  loop  module  to  SNPSAM. 

2.  Implement  other  changes  necessary  to  make  SNPSAM  fit  the  design  detailed 
in  (7:1-13). 

3.  Evaluate  the  performance  predicted  by  the  revised  version  of  SNPSAM. 

4.  Compare  the  performance  predicted  by  the  modified  program  with  the  perfor¬ 
mance  reported  by  General  Electric. 

While  the  above  were  the  formal  objectives  of  this  study,  the  following  objec¬ 
tives  give  it  practical  value  to  those  who  are  in  the  business  of  modeling  the  SP- 
100: 
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1.  Characterize  the  sensitivity  of  the  revised  program  to  changes  in  the  input 
file  and  in  the  source  code. 

2.  Demonstrate  methods  for  incorporating  changes  that  occur  in  the  SP-100. 

3.  Note  methods  that  facilitate  these  changes. 
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II.  Description  of  the  System  Model 


The  heat  rejection  system  of  the  SP-100  is  illustrated  in  Figure  2.1.  The 
salient  features  of  the  system  are  the  TEM  pump,  the  energy  conversion  assem¬ 
bly*,  the  radiator,  and  the  connective  piping.  Each  subsystem  is  modeled  in 
SNPSAM  and  will  be  described  in  the  following  sections.  First,  a  brief  description 
of  the  reactor  model  is  given. 


Figure  2.1.  Detail  of  the  Heat  Rejection  System  of  the  SP-100  (7:6) 


*  Unless  otherwise  stated,  "energy  conversion  assembly"  refers  to  the  primary  energy 
conversion  assembly,  not  the  energy  conversion  assembly  in  the  pump. 
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2.1.  Description  of  the  Reactor  Model 

The  reactor  model  consists  of  six  coupled  models:  a  kinetics  model,  a  reac¬ 
tivity  control  model,  a  reactivity  feedback  model,  a  decay  heat  model,  a  reactor 
thermal  model,  and  a  reactor  hydraulics  model.  All  the  models,  except  the  reactor 
hydraulics  model,  were  extracted  from  the  DSNP  library  module. 

The  kinetics  model  solves  the  neutron  point  kinetics  equations  using  six 
groups  of  delayed  neutrons.  The  reactivity  feedback  model  considers  the  reactivi¬ 
ties  induced  by  the  Doppler  effect,  core  expansion,  and  coolant  expansion.  The 
decay  heat  model  assumes  the  decay  heat  to  be  produced  only  by  gamma  radia¬ 
tion.  This  model  considers  up  to  five  groups  of  gamma  radiation.  The  core  ther¬ 
mal  model  computes  the  fuel,  cladding,  and  coolant  temperatures  within  the  reac¬ 
tor.  The  hydraulics  model  computes  the  pressure  loss  across  the  reactor. 

2.2.  Description  of  the  Primary  Heat  Exchanger  Model 

The  primary  heat  exchanger  is  a  flattened  section  of  conduit  designed  to 
effect  heat  transfer  from  the  fluid  flowing  through  it  to  a  bank  of  thermocouples. 
The  amount  of  heat  that  can  be  transferred  from  the  primary  heat  exchanger  to 
the  thermocouples  is  a  function  of  the  heat  transfer  coefficient  h  between  the  fluid 
and  the  heat  exchanger  material.  The  heat  transfer  coefficient  is  dependent  on 
fluid  properties  and  characteristics  of  the  flow,  both  of  which  are  dependent  on 
the  temperature  of  the  fluid. 

The  rate  of  heat  transfer  from  the  primary  fluid  to  the  thermocouples  is  also 
dependent  on  the  thermal- conductivity  k  of  the  intervening  materials  between  the 
fluid  and  the  thermocouples.  Thermal-conductivity  is  also  a  function  of  material 

temperature. 

Because  the  properties  of  the  materials  are  a  function  of  temperature,  a 
closed  solution  is  not  possible.  It  is,  therefore,  necessary  to  subdivide  the  fluid  in 


2.2 


the  heat  exchanger  into  discrete  segments.  The  temperature  throughout  each  seg¬ 
ment  is  a  constant.  This  allows  the  properties  of  each  segment  of  fluid  to  be 
determined  and  makes  it  possible  to  compute  the  heat  flow  out  of  the  segment. 


Figure  2.2  illustrates  how  the  coolant  flow  is  broken  into  segments.  The  tem¬ 
perature  at  the  inlet  to  the  heat  exchanger  is  T lf  the  second  segment  is  T2,  and 
so  on. 
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Figure  2.2.  Heat  Exchanger  Channel  Broken  into  Segments  (5:44) 

The  heat  transferred  into  or  out  of  an  arbitrary  segment  i  is 

Q  ~mcp  ( Ti+i—Tj )  (2.1) 

where  Q  is  the  heat  transferred  per  unit  time,  m  is  the  mass  flow  rate,  and  cp  is 
the  heat  capacity  of  the  coolant.  The  mass  flow  rate  is  assumed  to  be  constant 
throughout  the  coolant  loop  and  the  heat  capacity  is  a  function  of  the  coolant 
temperature  inside  the  segment.  Q  is  also  given  by 

Q=hA(Twaa-Taverage)  (2.2) 
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where 


T 

x  average 


Tj+i+Tj 

2 


(2.3) 


where  h  is  the  heat  transfer  coefficient  and  A  is  the  area  of  the  segment.  Com¬ 
bining  Eq  (2.1),  Eq  (2.2),  and  Eq  (2.3),  r,+1  may  be  written  as 

T,+l - l—\(mcp~^-)Ti+hAtmU) 

mc,+— 

The  old  and  new  heat  exchanger  designs  are  shown  in  Figure  2.3.  The 
configuration  of  the  two  designs  is  considerably  different.  The  former  design  was 
to  be  constructed  of  multiple  channels  of  varying  dimensions.  The  new  design  has 
only  one  channel  of  constant  dimension  (discounting  curvatures  at  the  inlets  and 
the  outlets). 


Figure  2.3.  Former  Heat  Exchanger  Design  (a),  and  the  New  (b)  (7:8  ;  14:57) 

The  heat  exchanger  assembly  is  symmetric  about  the  plane  that  halves  the 
primary  heat  exchanger.  It  is,  therefore,  possible  to  model  the  heat  exchanger  as 
one  cold  leg  and  one-half  the  hot  leg.  Because  the  heat  exchanger  is  modeled  with 
half  as  many  thermocouples  and  half  as  much  area  as  there  is  in  the  design,  the 
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computed  electrical  and  rejected  power  outputs  are  doubled. 

2.3.  Description  of  TEM  Pump  Model 

The  thermoelectric  electromagnetic  pump,  or  TEM  pump,  effects  a  pressure 
change  on  the  working  fluid  via  the  Lorente  force,  the  force  resulting  when  an 
electric  current  is  in  motion  through  a  magnetic  flux.  Because  the  electric  current 
is  in  motion  through  the  fluid,  the  force  acts  on  the  fluid.  The  electric  current  is 
provided  by  thermocouples  that  are  built  into  the  pump. 

The  force  acting  on  the  working  fluid  is  given  by 

F  -e(vXB)*  (2.4) 

The  configuration  of  the  TEM  Pump  is  illustrated  in  Figure  2.4.  The  B  field  is 
perpendicular  to  the  current  flow. 


Figure  2.4.  TEM  Pump  Configuration  (7:27) 


Eq  (2.4)  reduces  to 


*  Here,  bold  face  letters  represent  vectors. 
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F  =  a  I  B 


where  F  is  in  the  direction  of  fluid  flow  and  a  is  the  height  of  the  channel.  The 
pressure  generated  by  the  pump  is,  therefore, 

*f-TT-£T  (2-5) 

where  b  is  the  width  of  the  channel.  B  and  b  are  known.  The  total  current  /  is 
computed  by  summing  the  current  in  each  branch  of  the  equivalent  circuit  shown 
in  Figure  2.5 


Figure  2.5.  Equivalent  Circuit  for  the  TEM  Pump 


/*=/,+/„, +/y  (2.6) 

where  /  is  the  total  current,  Ie  is  the  current  passing  through  the  fluid,  Iw  is  the 
current  passing  through  the  wall  material,  and  Ij  is  the  current  passing  through 
materials  where  no  magnetic  field  exists.  The  electromotive  force  in  the  circuit  is 
given  by 

E-S-S.  (2.7) 
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where  Q  is  the  volumetric  flow  rate,  and  by 


E^IeRe-\-Ei=IwRw—lfRf  (2.8) 

where  Re,  Rw,  and  Rj ,  are  resistances  associated  with  the  coolant,  wall  material, 
and  absent  field  lines,  respectively.  Substituting  Eq  (2.7)  into  Eq  (2.8)  yields 

Iw—±-<I,R,+^-)  (2.9) 


and 

If  Eq  (2.9)  and  Eq  (2.10)  are  inserted  into  Eq  (2.6),  the  following  relation  is 
obtained 


I 


R  r  B  Q,K+R/  > 
RwRf  '  +  b  K  RwRf  ' 


(2.11) 


Solving  for  Ie 

RwRr  b  o  Rw+Rf 

/,«  — /■  / — ~r~^~(  (2.12 

R  b  R 

where  R  **Re  Rw  +RW  R j  +/? j  Re  . 

The  previous  TEM  pump  design  is  illustrated  in  Figure  2.6(a).  The  newest 
TEM  pump  design  is  illustrated  in  Figure  2.6(b).  Secondary  ducts  have  been 
added  such  that  the  thermoelectric  couples  are  fixed  between  the  primary  ducts 
and  the  secondary  ducts.  This  modification  effects  two  important  changes  in  the 
modeling  of  the  TEM  pump.  The  first  is  that  the  equivalent  circuit  for  the  pump 
is  altered,  the  second  is  that  the  boundary  condition  for  the  cold  interface  of  the 
thermoelectric  couples  is  no  longer  radiative,  but  is  a  constant  temperature  condi¬ 
tion. 


The  constant  temperature  condition  is  an  approximation.  Ideally,  the  condi- 
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tion  should  be  with  respect  to  heat  flow.  To  some  degree,  the  TEM  pump  acts 
like  a  heat  exchanger.  Since  the  primary  flow  is  in  close  proximity  with  (but  not 
counter-current  to)  the  secondary  flow,  an  amount  of  heat  will  be  transferred 
between  the  two  flows.  The  amount  of  heat  transferred  depends  on  the  length  of 
the  ducts  in  the  pump  and  the  heat  transfer  characteristics  of  the  fluid  and  pump 
materials.  In  this  study,  the  length  of  the  pump  was  assumed  short  enough  that 
the  heat  transferred  between  the  primary  and  secondary  is  negligible. 

The  equivalent  circuit  of  the  TEM  pump  is  modified  by  adding  circuit  ele¬ 
ments.  The  new  elements  are  Efecondary,  R“eonda' »,  R“condary,  and  Rf“°*d”y. 

When  these  elements  are  incorporated  into  the  equivalent  circuit,  the  pressure  rise 
across  the  TEM  pump  is  found  by  inserting  Eq  (2.12)  into  Eq  (2.5); 

\p_B  . _ B  Bw+Rf 

b  R  62  R  V 

where 

ft  ^ft  primary ft  secondary 


ft  x  ft  primary^. ft  secondary 


ftj  sft primary ft  secondary 


Q-Qprimary  ~^Qsecondary 


The  new  TEM  pump  configuration  is  symmetric  about  the  vertical  and  hor¬ 
izontal  center  lines  as  can  be  inferred  from  Figure  2.6(b).  It  is  possible,  therefore, 
to  model  the  TEM  pump  by  modeling  only  one  quadrant  of  the  pump. 

One  twelfth  of  the  mass  flow  rate  of  the  reactor  passes  through  each  half  of 
the  TEM  pump.  Each  primary  duct  is  symmetric  about  the  horizontal  center  line 
so  that  only  one  twenty-fourth  of  the  total  primary  mass  flow  rate  passes  through 
one  quarter  of  the  TEM  pump.  The  secondary  flow  is  split  by  the  two  separate 
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ducts  in  each  half  of  the  TEM  pump  so  that  one  twenty-fourth  of  the  total  secon¬ 
dary  mass  flow  rate  passes  through  one  quarter  of  the  TEM  pump. 

Since  the  primary  duct  is  identical  above  and  below  the  horizontal  line  of 
symmetry,  the  fluid  flow  in  both  halves  will  experience  an  equal  pressure  change. 
The  fluid  flow  in  the  upper  and  lower  secondary  ducts  will  also  experience  an 
equal  pressure  change.  In  the  adjacent  half  of  the  TEM  pump,  the  flow  is  return¬ 
ing  in  the  opposite  direction  such  that  the  total  pressure  change  induced  by  the 
entire  pump  is  equal  to  double  the  pressure  change  in  one  half  the  pump. 

The  new  boundary  condition  at  the  cold  interface  of  the  thermoelectric  cou¬ 
ples  is  implemented  by  setting  the  boundary  temperature  equal  to  the  secondary 
fluid  temperature.  The  radiative  condition  is  ignored  because  the  surface  area  is 
assumed  insufficient  to  radiate  enough  heat  away  to  cause  a  change  in  the  secon¬ 
dary  coolant  temperature. 

Figure  2.7  is  a  flow  diagram  for  the  heat  rejection  system.  It  indicates  one 
side  of  the  TEM  pump  handles  secondary  coolant  fluid  that  has  yet  to  enter  the 
radiator  while  the  other  side  handles  fluid  exiting  the  radiator.  This  suggests  that 
two  pump  modules  be  used,  one  that  uses  the  radiator  inlet  fluid  temperature  as 
the  boundary  condition,  and  a  second  using  the  radiator  outlet  fluid  temperature 
as  the  boundary  condition.  In  this  study,  however,  the  TEM  pump  is  modeled 
with  boundary  conditions  that  are  the  average  of  the  two  secondary  temperatures. 


2.4.  Description  of  the  Hydraulic  Model 

The  pressure  drop  in  a  pipe,  duct,  or  channel  is  given  by 

a p=-Lf£^ 

de  7  2  gc 


(2.13) 
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Figure  2.7.  Flow  Diagram  of  the  Heat  Rejection  System  (7:31) 

where  (/,=■--  (16:204),  w  and  t  being  the  width  and  thickness  of  the  channel, 
tu  I  t 

respectively.  /  is  the  friction  factor,  p  is  the  density  of  the  fluid,  v  is  the  velo¬ 
city  of  the  fluid,  and  ge  is  a  conversion  constant.  The  coolant  velocity  is 


m 

pwt 


Combining  Eq  (2.13)  and  Eq  (2.14)  yields  an  expression  for  A P. 


(2.14) 


2.5.  Description  of  the  Energy  Conversion  Model 

Thermoelectricity  is  generated  when  dissimilar  semiconductors  are  joined 
and  there  exists  a  temperature  gradient  perpendicular  to  the  joining  surface. 
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There  are  three  reversible  and  two  irreversible  phenomena  connected  with  the  gen¬ 
eration  of  thermoelectricity.  The  reversible  are  the  Seebeck  effect,  the  Thomson 
effect  and  the  Peltier  effect,  and  the  irreversible  processes  are  the  Fourier  effect  and 
Joule  heating  (17:155-156). 

The  amount  of  energy  converted  to  electricity  per  unit  temperature  difference 
is  given  by  the  Seebeck  coefficient 

or  **  lim  AQ/AT 
AT— *0 

where  A Q  is  the  change  in  heat  flow  and  AT  is  the  change  in  temperature. 

The  Peltier  effect  causes  heat  to  be  generated  at  the  junction  of  different 
materials  that  are  at  the  same  temperature  when  an  electrical  current  is  present. 
The  heat  generated  is  proportional  to  the  Peltier  coefficient  npn ,  given  by 

Qp^pn1 

where  Q  is  the  heat  flow  across  the  junction,  I  is  the  electrical  current,  and  where 

Kpn  “  np  -*n  “  np  +7rn 

where  np  and  irn  are  the  Peltier  coefficients  for  the  p  and  n  materials,  respec¬ 
tively. 

The  Thomson  effect  causes  heat  to  be  generated  in  a  homogeneous  material 
where  there  exists  a  temperature  gradient  when  an  electric  current  is  present.  The 
definition  of  Thomson  coefficient,  r,  is 

r*  lim 
AT— *0 

Joule  heating  occurs  when  an  electrical  current  passes  through  a  material. 

The  heat  evolved  is  equal  to 

Q,  -  PR 


A  OR 
A  T 
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where  R  is  the  resistance  presented  by  the  material. 


The  diffusion  of  heat  in  a  solid  that  is  experiencing  a  temperature  gradient  is 
described  by  the  Fourier  effect  where  the  heat  flow  is  given  by 

,dT 
q  *=  ~k— 
dx 

where  q  is  the  heat  flux  and  k  is  the  thermal-conductivity. 

The  coefficients  of  the  three  reversible  process  are  related  by  Kelvin’s  rela¬ 
tions  (4:4).  Kelvin’s  relations  are 

7T  -  ocT  (2.15) 

and 

,da 


r«  T- 


Combining  Eq  (2.15)  and  Eq  (2.16)  yields 


dT 


(2.16) 


dT  dT  T 


The  energy  balance  as  a  function  of  time  within  the  thermoelectric  material  is 
given  by 


S  lSt  .  .  nu\&T  r,f>T 

+«(0  +  G(0-5-  -<tt 


(2.17) 


where 


Q{t)  =  pJ\t) 
G{t)=Tj(t) 


and 
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where  w  is  the  coolant  density,  satisfying  the  boundary  conditions 

mo-uo.  T(L,t)-TcM(t) 

and  the  initial  condition 

T(x, 0)  *  T0(x) 

All  terms  of  Eq  (2.17)  are  moved  to  the  left  hand  side.  Then,  each  side  of 
the  resulting  equation  is  multiplied  by  a  test  function,  v,  and  integrating  over 
space. 


V^xkl!z  + 


ST  nST  1. 

- - t;C—  MX  *=  0 

6x  St  f 


(2.18) 


Integrating  the  first  term  of  Eq  (2.18)  by  parts  over  the  interval  from  0  to  L 

rL  S  t,ST v,  /,  £7\  | £,  ST  ,  fn  . 

I  ”&<*■&)*  - *■(*-&) l»  - 1  kteJTix  (2'1J> 

The  first  term  on  the  right  hand  side  of  Eq  (2.10)  is  zero  because  v  is,  by 
definition,  equal  to  zero  at  0  and  L.  Substituting  the  second  term  in  Eq  (2.19) 
into  Eq  (2.18)  and  multiplying  both  sides  of  the  resulting  equation  by  —1 


J t 


ST  nST\, 

~Sx  vC~sT  j  ™  0 


(2.20) 


The  Galerkin  approximation  suggests  the  solution  of  (2.20)  has  the  form 


N 


T-ETjitftix) 

y- 1 


AT 


v "  £«*(«)*.•(*) 


(2.21) 


(2.22) 


i-l 


where  <t>i(x),  i**l,2,...,yV ,  are  basis  functions  (8:190-192).  Eq  (2.21)  and  Eq  (2.22) 
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are  substituted  into  Eq  (2.20)  to  obtain 
N  L 


N 

Eif(0 

«-i 


+  p^L k±4,,(I)±4,j(I)Ti(t)dx  - 


N  _L 


l  Q(t)<t>i{x)dx  -  G ( t )<f>i {x)-j-<j)j{x) Tj (t)dx 


(2.23) 


where  Tj(t)  —  —Tj(x).  Eq  (2.23)  is  the  form  of  Eq  (2.20)  to  be  solved  numeri- 


N 


cally.  Matters  are  simplified  if  Eq  (2.23)  is  multiplied  through  by  t \(t)  and 

i-i 


each  term  studied  separately.  The  first  term  is 

N  N  L 


X>.(0El  ^.(*)<Mx)7V(0<fr 

.-1  y-1 


(2.24) 


Carrying  out  the  summations  in  the  first  first  term 


+  lLc^(x^x)f^t)dx  +  ...  +  lLcu*)M*)TN(  t)<fc]  + 
”2(0  [£c<U*)*x(x)Tx(t)dx  +  £c+Ax)Hx)T4t)dx  +  ...  +  JtiC,^*)M*)TJV(t)dx]  + 

®w(o[X  C'^/v(*)^t(*)T,(0‘k  +  X  C^N{x)d2(x)T^t)dx  +  ...  +  £  C^(z)^(i)7V(<)rfx] 
The  above  summation  can  be  written  as  a  matrix  product 

«,(*) 

va(x) 

*>ti(x) 

The  Vf(x)  vector  appears  in  each  term  of  the  matrix  form  of  Eq  (2.23).  Since  the 


X  Cd>i(x)<l>i{x)dx 


X  Cdy(x)d,(x)dx  . 


X  C0J (x)<f>N(x)dx 


X  C<t>N(x)dN(x)dx 


r,(<) 

W 

rN(t) 


equation  is  set  equal  to  zero,  the  vi  (x )  vector  cancels. 

The  basis  functions  are  chosen  to  be  piecewise  linear  shape  functions,  or  "hat" 


functions,  as  in  Figure  2.8  (11:472-476).  These  functions  are  defined  as 

t 

l/h(x-xn)+l  xe[xn_vxn] 

<M*)  =  ■  -1/A(x-Xn)+1  x  e  [xB,xB+1] 

0  x  <  [x„_i,xB+1] 

where  h  defines  the  slope  of  the  hat  function. 


y 

y 

j- 2  j- 1  J- 

ay  j+ 1 

Figure  2.8.  "Hat"  Functions  Used  in  Solving  Eq  (2.23)  (11:473) 

Each  element  of  the  matrix  can  now  be  determined.  Elements  possessing  the  form 
/  N  C<t>„dx  have  the  solution 


/  NC<ftdx  *  J  ’  4-(x— xn)+l  dx 

*0  J  (l 


Letting  u  =  -i-(x— x„)+l 
h 


cC ~  ~  t‘ 


Elements  possessing  the  form  J  <t>n-\ have  the  solution 

Xq 


is"  C<t>n_x<t>ndx  **  jt^C<t>n-\<i>ndx 
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1 


jr(*-*n)+l 

Again,  letting  o  ■»  ^-(i  —xn  )+l 
n 

J  " C<f>n_x<f>ndx  —  Ch  f  (1— u)udu  *=  —Ch 

*o  6 

All  other  elements  are  zero.  The  matrix  is,  therefore,  a  tridiagonal  matrix. 

Each  element  may  be  written  more  compactly  as 

cij  =  l  C 4>i(x)<f>j(x)dx 

Each  of  the  other  terms  in  Eq  (2.23)  may  also  be  written  in  this  form 

Qi  =  l  Q(t)<t>i(x)dx 
Gii  s  )-^Mx)dx 

Each  of  these  terms  may  be  solved  for  by  using  the  same  procedure  outlined  for 
the  first  term. 

In  matrix  notation*,  Eq  (2.23)  becomes 

[C] T  +  \K) T  -  Q  -  (G]T  =  0  (2.25) 

Letting  \K ]  —  (G]  *=  [//]  and  moving  Q  to  the  right  hand  side,  Eq  (2.25)  becomes 

[C]T  +  [H] T  *  Q  (2.26) 


f'C 

Jx«-1 


*  Capitol  letters  enclosed  by  brackets  are  used  here  to  represent  a  matrix.  Bold  face  letters 
represent  column  vectors. 
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Using  a  finite  difference  scheme,  T  is  written  as  the  backward  difference 


Eq  (2.26)  may  then  be  written  as 

[[C]  +  A*{/f]  Jr*+1  =  A«Q*+1  +  [C]T*  (2.27) 

All  factors  of  Eq  (2.27)  are  known  except  T*+1.  A  suitable  matrix  algorithm  can 
solve  for  Ti+1,  as  done  in  SNPSAM. 

2.6.  Description  of  the  Secondary  Heat  Exchanger  Model 

On  the  opposite  side  of  each  bank  of  thermocouples  is  the  secondary  heat 
exchanger,  making  up  the  cold  leg  of  the  energy  conversion  assembly.  It  is  identi¬ 
cal  in  every  way  to  the  primary  heat  exchanger  except  that  fluid  flow  is  counter 
current  to  the  primary  flow. 

The  input  parameter  to  the  secondary  model  is  the  heat  flow  Q  from  the 
thermocouple  model.  Therefore,  the  temperature  of  the  fluid  in  each  segment  is 

T,+I  -  ~S~+Ti 

mcp 

2.7.  Description  of  the  Radiator  Model 

Coolant  fluid  enters  the  radiator  channel  at  the  secondary  heat  exchanger 
exit  temperature  and  leaves  the  radiator  channel  at  the  secondary  heat  exchanger 
inlet  temperature.  The  method  for  computing  the  temperature  along  the  radiator 
channel  is  the  same  method  used  to  compute  the  temperature  along  the  primary 
heat  exchanger,  except  that  the  wall  temperature  is  computed  based  on  the  heat 
flow  permitted  by  the  radiator. 

The  radiator  is  modeled  as  a  black-body  emitter  where  the  heat  leaving  from 
its  surface  is  governed  by 
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"1 


F-tj  mr'-T*') 

where  F  is  the  heat  flux,  e  is  the  emissivity  of  the  radiator,  and  /  is  the  view  fac¬ 
tor  for  the  radiator.  T  in  the  above  equation  is  Tcold(t),  the  boundary  condition 
at  the  cold  junction  of  the  thermocouples  (as  described  in  section  2.5).  e  and  / 
are  defined  in  the  input  file  SYS.INP  under  RADIATOR  PARAMETERS  (see 
Appendix  A). 


I 

i 
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HI.  Procedure 


Two  hundred  and  seventy  lines  of  Fortran-77  instruction  were  added  to 
SNPSAM.FOR.  The  revision  has  been  given  the  version  name  Rev  GNE88M-X. 
Eighty  eight  lines  have  been  added  to  SYS.INP.  The  revision  of  SYS.INP  is 
reproduced  in  this  report.  Only  the  major  subroutines  that  were  added  to 
SNPSAM  have  been  appended  to  this  document.  Some  of  the  changes  are  dis¬ 
cussed  in  the  following  sections.  Other  changes,  such  as  COMMON,  WRITE,  and 
READ  statements,  for  example,  were  added  to  SNPSAM  primarily  to  handle 
changes  in  data  structure  and  will  not  be  discussed. 

SNPSAM  is  sensitive  to  changes,  whether  these  changes  are  made  in  the 
source  code,  or  in  the  input  parameters.  Specifically,  the  thermocouple  model, 
TMSTR,  uses  an  iterative  method  that  fails  to  converge  for  input  that  is  outside 
a  narrow  range*.  To  accommodate  this  problem,  a  procedure  was  used  where 
changes  were  made  incrementally.  The  drawback  of  this  method  is  that  it  requires 
extra  steps  and  is  time-consuming.  A  method  of  predicting  if  a  given  change  will 
cause  the  execution  of  SNPSAM  to  fail  was  not  developed. 

The  procedure  used  to  change  input  parameters  is  illustrated  by  an  account 
of  how  the  heat  exchanger  configuration  was  modified  in  the  following  section. 

The  procedure  used  to  make  changes  to  the  source  code  is  discussed  in  section 
3.1.2. 

3.1.  Making  Changes  to  Input  Parameters 

Input  parameters  are  changed  by  editing  the  input  file  SYS.INP.  SYS.INP  is 
a  formatted  input  file  that  is  self-explanatory.  A  brief  description  of  each  input 

*  The  solution  space  was  not  determined  in  this  study. 
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parameter  is  given  next  to  the  field  where  its  value  is  entered.  An  example  input 
file  (SYS.INP  Rev  GNE88M-X)  is  given  in  Appendix  A.  Pipe  and  TEM  pump 
dimensions  are  examples  of  input  parameters  that  were  modified  but  did  not 
require  incremental  changes.  That  is,  SNPSAM  successfully  executed  when  these 
changes  were  made.  Other  parameters  were  not  amenable  to  this  straight-forward 
approach.  Heat  exchanger  dimensions  and  radiator  surface  area  are  examples  of 
parameters  that  had  to  be  changed  incrementally.  The  number  of  increments 
required  varied  depending  on  the  parameter.  Ten  increments  was  typically  neces¬ 
sary. 

The  parameters  used  for  SNPSAM  to  simulate  the  particular  SP-100  design 
were  fund  in  (7:1-13),  but  a  few  elements  of  this  data  are  suspect.  These  data 
were  "eyeballed"  from  illustrations  obtained  on  the  design.  No  numerical  data 
were  available  in  these  cases.  The  "eyeballed”  parameters  include  the  length  of  the 
radiator  coolant  channel,  some  dimensions  of  the  TEM  pump,  and  the  intercon¬ 
nective  piping  into  and  out  of  the  TEM  pump  and  heat  exchangers.  The  error 
associated  with  the  radiator  length  and  the  TEM  pump  dimensions  is  approxi¬ 
mately  two  and  one  percent,  respectively.  The  error  in  the  interconnective  piping 
lengths  was  not  determined*. 

3.1.1.  Changing  the  Heat  Exchanger  Geometry.  The  new  heat 
exchanger  was  simulated  as  a  longer  one-channel  version  of  the  old  heat 
exchanger.  The  input  file  made  available  for  this  study  used  only  one  segment  for 
each  channel.  The  initial  attempt  to  simulate  the  new  heat  exchanger  using  only 
one  segment  failed.  It  was  necessary  to  use  more  than  one  segment  for  the  one 
channel  heat  exchanger  model  (as  described  in  section  2.2). 

Given  an  input  file  that  simulates  n  segments,  the  procedure  to  generate  an 

*  No  accurate  scale  for  the  lengths  was  found  and,  furthermore,  the  figure  was  rotated  on 
two  axis. 
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input  file  that  models  n+1  segments  requires  that  one  set  of  initial  conditions  be 

added  to  the  input  file.  The  steps  necessary  to  compute  a  new  set  of  initial  condi¬ 
tions  are  as  follow: 

1.  Indicate  in  the  input  file  that  one  more  segment  is  going  to  be  used. 

2.  Append  a  new  line  of  initial  conditions  (i.e.,  a  new  set  of  initial  conditions  for 
another  segment)  by  duplicating  the  last  line  of  initial  conditions. 

3.  Run  SNPSAM  in  the  steady-state  mode  until  a  steady-state  solution  is 
found.  The  input  file  to  SNPSAM,  SYS.INP,  offers  an  option  of  steady-state 
operation  or  transient  operation.  In  the  steady-state  mode,  SNPSAM  simu¬ 
lates  the  operation  of  the  SP-100  at  a  constant  power  level.  If  the  initial  con¬ 
ditions  are  consistent  with  the  steady-state  operation  at  that  power,  no 
change  in  any  of  the  outputs  will  be  observed  in  that  particular  steady-state 
simulation.  However,  if  the  initial  conditions  are  inconsistent  with  the 
steady-state  operation  at  that  power,  the  outputs  of  the  simulation  will 
change  until  they  converge  to  steady-state  levels. 

4.  Call  the  output  file  that  contains  the  new  initial  conditions  and  append  it  to 
the  input  file.  SNPSAM  writes  the  steady-state  initial  conditions  to  two  files 
after  it  has  been  executed.  These  files  are  named  SSTX.INP  and  SSTP.INP. 
SSTX.INP  is  the  output  file  containing  initial  conditions  for  the  primary 
energy  conversion  assembly.  SSTP.INP  is  the  output  file  containing  initial 
conditions  for  the  pump  energy  conversion  assembly. 

5.  Edit  the  input  file  by  removing  the  old  initial  conditions  and  replace  them 
with  the  appended  initial  condition. 

6.  Repeat  steps  one  through  five  until  the  needed  number  of  initial  conditions 
exists  (i.e.,  the  number  of  segments  used  in  the  simulation). 

This  procedure  was  repeated  until  ten  segments  were  represented  in  the  input 
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3.1.2.  Changing  Other  Input  Parameters.  Fortunately,  the  procedure 
required  to  change  most  input  parameters  is  not  as  complicated  as  the  one 
described  above.  Other  parameters  do  not  require  additional  segments  be  com¬ 
puted  and  entered  into  the  input  file.  To  increase  the  reactor  power  from  2000 
kW  to  6800  kW  as  done  in  this  study,  for  example,  it  was  necessary  only  to  com¬ 
pute  one  intermediate  set  of  initial  conditions.  When  6800  was  entered  as  the 
reactor  power,  having  made  no  other  changes,  SNPSAM  failed  upon  execution 
(due  to  the  sensitivity  of  the  initial  conditions).  The  intermediate  step  was  to 
enter  4000  as  the  reactor  power,  execute  SNPSAM  in  steady-state  mode,  retrieve 
the  initial  conditions  written  into  SSTX.INP  and  SSTP.INP,  enter  these  condi¬ 
tions  into  SYS.INP  in  place  of  the  former  conditions,  and,  finally,  execute 
SNPSAM  with  6800  entered  in  SYS.INP  as  the  reactor  power. 

3.2.  Modification  to  the  SNPSAM  Source  Code 

Entering  changes  into  the  source  code  may  have  the  same  effect  as  changing 
initial  conditions  in  the  input  file  (obviously  a  WRITE  statement  and  similar 
instructions  will  not  have  this  result).  For  example,  in  this  study  it  was  necessary 
to  reduce  the  mass  flow  rate  in  the  secondary  and  primary  heat  exchangers  by 
half.  To  do  so,  instead  of  dividing  by  2,  the  mass  flow  rate  was  divided  by  a  tem¬ 
porary  variable.  The  variable  was  included  in  the  input  file  and  SNPSAM  was 
modified  so  that  it  would  be  read.  The  variable  was  initially  given  a  value  of  one 
and  then  was  varied  incrementally  as  described  in  the  above  section.  The  major 
changes  made  to  the  source  code  are  discussed  below. 

3.2.1.  Addition  of  the  Secondary  Model.  A  subroutine  that  simulates 
the  operation  of  the  secondary  was  written,  tested,  and  appended  to  SNPSAM. 
The  subroutine  was  titled  SECONDARY.  A  subroutine  that  simulates  the  opera¬ 
tion  of  the  radiator  was  written,  tested,  and  appended  to  SNPSAM.  This 


subroutine  was  titled  RADIATOR.  SECONDARY  and  RADIATOR  are  coupled 
subroutines  that  compute  the  boundary  conditions  for  the  thermocouple  model. 

In  Appendices  C  and  D  is  a  listing  of  SECONDARY  and  RADIATOR,  respec¬ 
tively. 

3.2.2.  Modifications  to  the  TEM  Pump  Model.  Since  the  TEM  pump 
operates  on  the  primary  and  secondary  (i.e.,  the  primary  and  secondary  are  cou¬ 
pled  to  thermocouples  that  drive  the  pump),  a  separate  subroutine  accounting  for 
the  secondary  channels  in  the  pump  could  not  be  written.  Therefore,  equations 
within  the  TEM  pump  subroutine,  titled  EMPMP,  were  modified  to  account  for 
secondary  channels  in  the  pump.  The  equations  are  embedded  in  a  listing  of 
EMPMP,  which  has  been  reproduced  in  Appendix  E.  The  nature  of  these  changes 
are  discussed  in  section  2.3. 

3.2.3  Addition  of  Support  Subroutines.  The  subroutine  responsible  for 
computing  the  Nusselt  number,  CONVEC,  was  duplicated.  The  duplicate  was 
modified  and  entered  into  SNPSAM  following  the  unmodified  subroutine.  This 
was  necessary  because  the  choice  of  COMMON  variables  in  SNPSAM  does  not 
allow  CONVEC  to  be  called  by  both  the  primary  heat  exchanger  model  and  the 
secondary  heat  exchanger  model.  The  added  subroutine  is  titled  CONVECSEC. 
CONVECSEC  is  identical  to  CONVEC  except  that  COMMON  statements  con¬ 
taining  information  about  the  primary  were  changed  to  contain  information  about 
the  secondary. 

The  above  procedure  was  also  performed  to  the  subroutines  DPIPE  and 
HYDR.  DPIPE  computes  the  pressure  loss  in  the  primary  piping  while  HYDR 
computes  the  mass  flow  rate  in  the  primary.  After  changes  were  made  to  copies 
of  these  subroutines,  they  were  resubmitted  under  the  names  DPIPESEC  and 
HYDRSEC.  These  subroutines  are  called  by  SECONDARY. 


3.3.  Role  of  Secondary  in  Coupled  Model 

The  reactor  model  is  given  the  value  of  TECI,  the  reactor  coolant  tempera¬ 
ture  at  the  reactor  inlet,  and  returns  the  value  of  TECX,  the  reactor  coolant  tem¬ 
perature  at  the  reactor  exit.  The  first  value  of  TECI  that  is  passed  to  the  reactor 
model  is  chosen  by  the  user  as  an  initial  condition.  At  all  other  times,  TECI  is 
computed  by  the  primary  heat  exchanger  module,  HTEX.  HTEX  generates  a 
value  for  TECI  by  using  TECX  as  input. 

The  heat  exchanger  module  modifies  TECX  in  geometric  intervals  along  the 
axial  length  of  the  heat  exchanger.  These  intervals  are  termed  segments.  While 
at  a  segment  *  of  known  temperature  T,  ,  the  temperature  of  the  following  seg¬ 
ment,  T,+1,  is  estimated.  The  heat  flow,  Q,  out  of  this  segment  is  given  by 
<j>=mcp(7,J+1--7,l).  Q  for  each  segment  is  needed  by  TMSTR,  the  module  tasked 
to  simulate  the  thermoelectric  converters.  TMSTR  returns  a  value  for  the  tem¬ 
perature  of  the  primary  heat  exchanger  and  thermoelectric  interface,  Thot .  HTEX 
uses  Tkot  to  compute  the  value  for  Tl+1. 

To  compute  T*,  TMSTR  must  guess  the  temperature  at  the  interface  of  the 
thermoelectric  converters  and  the  secondary  heat  exchanger.  The  guess  is  termed 
T*old .  The  secondary  module,  SECONDARY,  modifies  the  guess  to  account  for 
how  much  heat  that  can  be  rejected  from  the  radiator.  The  modified  guess  is 
termed  T?0fdl.  TMSTR  then  recomputes  Tkot  and  T?old  until  |  T?0fd—  T?0 ld  |  is 
less  than  e,  some  convergence  criteria. 

The  heat  flow  across  the  thermoelectric  converter  and  secondary  heat 
exchanger  interface  is  required  by  SECONDARY  to  compute  Tc"  J1.  SECON¬ 
DARY  begins  at  the  cold  leg  inlet  with  an  initial  guess  for  the  inlet  temperature, 
TESI* .  SECONDARY  advances  segment  by  segment  down  the  secondary  chan¬ 
nel  computing  the  average  temperature  of  the  working  fluid  in  that  segment  based 
on  the  heat  flow  into  the  channel  from  the  ther>noelectrics.  The  last  segment  is  at 


the  outlet  of  the  secondary  heat  exchanger.  The  fluid  enters  the  radiator  at  the 
secondary  exit  temperature,  TESX.  SECONDARY  also  calls  HYDRSEC,  the 
hydraulic  model  for  the  secondary  coolant  loop,  which  returns  a  value  for  the 
mass  flow  rate  within  the  secondary. 

The  module  RADIATOR  accepts  TESX  as  the  fluid  temperature  of  the  first 
segment  in  the  radiator  channel  and  computes  successive  fluid  temperatures  of 
each  segment.  The  temperature  of  successive  segments  decreases  because  there  is 
heat  loss  out  of  the  segment  by  radiative  heat  transfer. 

The  outlet  temperature  of  the  radiator  is  TESIn+1.  This  becomes  the  new 
secondary  inlet  temperature  and  is  compared  with  the  old  inlet  temperature, 
TESIn .  If  they  are  within  the  convergence  criteria,  then  TMSTR  is  passed  the 
boundary  value, 
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IV.  Results 


4.1.  Temperature  and  Heat  Flow  Profiles  of  the  Secondary 

The  temperature  profiles  exhibited  by  the  secondary  are  shown  in  Figure  4.1. 
The  profiles  were  extracted  from  a  steady-state  run  of  SNPSAM  where  the  reactor 
thermal  power  was  set  to  6.8  megawatts,  but  the  profiles  are  typical  of  those 
observed  in  this  study.  The  heat  flow  per  segment  from  the  primary  to  the  secon¬ 
dary  is  shown  in  Figure  4.2  as  a  function  of  relative  length  along  the  heat 
exchanger. 

4.2.  Energy  Balances 

4.2.1.  Steady-State  Energy  Balance.  Figure  4.3  reports  the  various  heat 
flow  rates  that  occur  in  the  SP-100.  The  reactor  thermal  power,  the  thermal 
power  into  the  energy  conversion  assembly,  the  electrical  power  generated  by  the 
energy  conversion  assembly,  and  the  thermal  power  dissipated  by  the  radiative 
panel,  is  presented.  The  purpose  of  the  figure  is  to  provide  an  energy  balance  of 
the  system.  The  mode  of  operation  is  steady-state.  In  theory,  the  power  rejected 
added  to  the  electric  power  should  equal  the  power  input  to  the  thermocouples. 
The  difference  is  as  much  as  700  kilowatts,  representing  approximately  eleven  per¬ 
cent  of  the  total  power. 

It  is  evident  from  Figure  4.3  that  the  initial  conditions  for  the  secondary  are 
not  consistent  with  steady-state  operation  at  that  reactor  thermal  power,  indi¬ 
cated  by  the  varying  power  level,  but  the  rejected  thermal  power  settles  to 
steady-state  values  in  approximately  ten  seconds. 
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Function  of  Relative  Length 


Power  (kilowatts) 


4.2.2.  Transient  Energy  Balance.  Figure  4.4  reports  thermal  powers  in 
the  SP-100  during  a  50  second  simulation  where  the  reactivity  is  linearly  ramped 
from  0.0  to  1.0a:  10-4  (absolute  reactivity),  as  in  Figure  4.5.  The  choice  of  the 
ramp  function  and  its  magnitude  was  purely  arbitrary.  Similar  results  are 
obtained  for  other  reactivity  insertions.  An  example  of  a  negative  reactivity  inser¬ 
tion  is  included  in  Appendix  G.  The  power  balance,  reactivity  insertion,  tempera¬ 
tures,  and  mass  flow  rates  are  represented  in  Figures  G.l,  G.2.,  G.3,  and  G.4, 
respectively. 

Temperatures  at  key  points  in  the  system  as  they  vary  with  time  are  shown 
in  Figure  4.6.  The  mass  flow  rates  of  the  primary  and  secondary  as  a  function  of 
time  are  shown  in  Figure  4.7. 

The  rejected  power  is  seen  to  oscillate  in  the  first  few  seconds  of  simulation, 
after  which  it  becomes  smoothly  varying.  The  difference  between  the  sum  of  the 
electrical  power  and  the  rejected  power  with  the  power  entering  the  thermocouples 
is  approximately  200  kilowatts,  or  about  five  percent  of  the  total  thermal  power 
(not  taking  into  account  the  differences  that  occurs  where  the  oscillations  exist). 

Figure  4.8  reports  the  power  balance  results  of  a  quasi-steady-state  simula¬ 
tion  using  the  transient  option.  To  accomplish  this,  5xl0-6  absolute  reactivity 
was  inserted  in  a  linear  ramp  function  shown  in  Figure  4.5.  The  power  was  then 
allowed  to  reach  steady-state  values.  The  power  leveled  out  at  6003.6  kW, 
approximately  1.50%  greater  than  the  total  reactor  power. 

The  rejected  plus  electrical  power,  indicated  in  Figure  4.8,  is  closer  in  value 
to  the  input  power  to  the  thermocouples  than  the  steady-state  solution,  even 
though  they  operated  at  close  to  the  same  reactor  thermal  power. 


4.5 


Power  (kilowatts) 


React  (absolute) 


Temp  (K) 


Mass  Flow  (Kg/sec 


8 

oi 


8 

ty 


8 

© 

6 


- ! - [ 

- , - , - 

Secondary 

_ 1 - ! - 

Primary 

_ 1 _ 1 _ 

0.000 


tooo  aooo  3.000 

Time  (seconds) 


4.000 


xlO 


5.000 


Figure  4.7.  Mass  Flow  Rate  in  the 
tion  of  Time  for  Transient  Case 


Primary  and  Secondary  as  a  Func- 
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4.3.  Comparison  with  General  Electric  Results 

Results  obtained  from  SNPSAM  (Rev  GNE88M-X)  and  results  reported  by 
General  Electric  are  presented  in  Table  4.1.  The  results  were  generated  for  a  reac¬ 
tor  thermal  power  of  6.8  megawatts,  steady  state  operation.  The  numbers  quoted 
for  SNPSAM  are  taken  from  the  quasi-steady-state  case  of  Figure  4.8  immediately 
following  the  oscillations. 

The  overall  difference  between  the  two  is  that  SNPSAM  predicts  consistently 
lower  temperatures  and  higher  electrical  power  output  than  predicted  by  General 
Electric.  The  mass  flow  rate  generated  by  SNPSAM  is  also  lower  than  reported 
by  Genera!  Electric. 


Table  3.1.  SNPSAM  Results  Compared  with  General  Electric  Results 


Parameter 

SNPSAM 

General  Electric 

Reactor  Thermal  Power 

6.8  MW 

6.8  MW 

Power  Rejected 

6.28  MW 

6.-15  MW 

Radiator  Temperature 

76-4  K 

774  -  812  K 

Mass  Flow  Rate  Per  Module 

1.86kg/see 

2.12  kg/sec 

Loop  Pressure  Drop 

17  kPa 

24  -  52  kPa 

Electrical  Power 

■119  kW 

336  kW 

Voltage 

170  V 

200  V 

Fluid  Temperature 

1288  K 

1320  K 

Hot  Junction  Temperature 

12-19  K 

1275  K 

Cold  Junction  Temperature 

791  K 

855  Iv 

Conversion  Efficiency 

6.18 
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V.  Conclusions  and  Recommendations 

This  study  succeeded  in  modeling  a  secondary  coolant  loop  in  the  SP-100. 
The  accuracy  of  the  model  is  evidenced  by  the  following: 

■  energy  balance  of  the  system 

■  temperature  profiles  in  the  system  are  characteristic  of  like  systems 

■  results  compare  reasonably  well  with  results  reported  by  General  Electric 

Because  the  error  in  the  steady-state  simulation  mode  was  not  found,  the  results 
used  to  base  the  success  of  the  secondary  model  are  those  of  the  quasi-steady-state 
case.  It  is  reasonable  to  believe  that  if  the  steady-state  simulation  mode  func¬ 
tioned  correctly,  it  would  report  results  that  match  the  results  reported  by  the 
quasi-steady-state  case  (i.e.,  slowly  varying  transient  case). 

Figures  4.4,  4.8,  and  G.l  indicate  that  an  energy  balance  exists  in  the  system. 
In  each  case,  the  flow  into  the  thermocouples  is  matched  by  the  sum  of  the  electri¬ 
cal  power  and  the  rejected  power. 

Table  4.1  indicates  reasonable  agreement  between  SNPSAM  (Rev  GNE88M- 
X)  and  General  Electric.  Notable  for  close  agreement  are  the  rejected  power,  the 
radiator  temperature,  the  primary  fluid  temperature  (average),  and  the  hot  junc¬ 
tion  temperature.  These  output  parameters  are  within  3%*  of  values  reported  by 
General  Electric.  The  mass  flow  rate  per  module,  the  loop  pressure  drop,  the 
electrical  power,  the  voltage,  and  the  conversion  efficiency  reported  in  Table  4.1 
represent  12,  29-67,  25,  15,  and  14  percent  disagreements,  respectively,  with  Gen¬ 
eral  Electric. 

*  Relative  error. 
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The  temperature  and  heat  flow  profiles  generated  by  the  secondary  model  are 
consistent  with  typical  counter-current  heat  exchangers  (16:406-411).  An  almost 
constant  temperature  difference  between  the  hot  leg  and  the  cold  leg  is  established 
(Figure  4.1).  The  heat  flow  rate  profile  associated  with  the  aforementioned  tem¬ 
perature  profile  (Figure  4.2)  is  also  nearly  constant. 

The  source  of  the  power  imbalance  that  occurs  in  the  steady-state  simulation 
(Figure  4.3)  was  not  pinpointed.  This  imbalance  occurred  to  a  lesser  degree  in  the 
unrevised  version  of  SNPSAM,  so  it  is  safe  to  assume  the  problem  is  rooted  in  the 
subroutines  that  were  added  or  in  the  interaction  between  the  subroutines  that 
were  added  and  the  subroutines  that  are  in  the  unrevised  SNPSAM. 

The  source  of  the  power  imbalance  is  likely  to  be  the  interaction  between  the 
thermoelectric  model  and  the  secondary  model.  Why  the  imbalance  occurs  is  attri¬ 
buted  to  the  inability  of  the  thermocouple  model  to  compute  a  solution  based  on 
the  information  given  to  it  by  the  secondary  model.  This  is  evidenced  by  Figure 
4.3  and  Figure  4.8.  When  these  figures  are  superimposed,  the  first  few  seconds  of 
the  rejected  power  for  each  case  match  identically.  After  the  first  few  seconds, 
however,  the  two  curves  separate.  It  appears  the  quasi-steady-state  curve  recovers 
from  an  oscillation  while  the  steady-state  case  incorrectly  converges  to  a  lower 
value.  Since  the  finite  element  method  is  used  in  the  thermocouple  model,  a  possi¬ 
ble  remedy  for  the  problem  is  to  replace  the  piecewise  linear  shape  functions  with 
piecewise  quadrature  functions.  Also,  higher  order  terms  could  be  added  to  the 
finite  difference  that  defines  T  (section  2.5). 

An  idiosyncrasy  developed  having  to  do  with  generating  an  output  file. 

When  75  and  200  seconds  were  used  as  simulation  times,  PLOT.TMP  was  opened 
but  no  output  data  was  written  to  the  file.  Other  simulation  times,  20,  50,  125, 
and  250  seconds,  for  example,  did  not  result  in  this  error. 
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The  subroutine  that  writes  PLOT.TMP,  PLOTE,  was  modified  in  this  study  so  it 
is  likely  that  this  is  the  source  of  damage. 

The  uncertainty  associated  with  the  interconnective  piping  lengths  (section 
3.1.)  were  not  a  source  of  significant  error.  For  example,  when  the  lengths  of  pip¬ 
ing  connecting  the  heat  exchanger  with  the  TEM  pump  are  doubled  (i.e.,  assume 
100%  error  in  the  known  length),  the  corresponding  change  in  A P  is  negligible  ( 
<  2%  ).  This  is  true  because  the  lengths  of  piping  connecting  the  heat 
exchangers  are  much  shorter  than  those  transporting  coolant  to  the  radiator 
panels. 

General  Electric  considers  the  heat  loss  caused  by  conduction  between  the 
channels  carrying  fluid  to  the  radiator  and  channels  carrying  fluid  returning  from 
the  radiator.  The  channels  are  counter- current  to  each  other  and  are  in  close 
proximity  (see  Figure  2.7),  causing  heat  to  be  exchanged  efficiently  between  the 
two  flows.  Consequently,  the  efficiency  of  the  radiator  is  reduced.  Also,  General 
Electric  considers  the  heat  lost  from  small  holes  in  these  channels  used  for  thaw¬ 
ing  the  coolant  during  start-up  (3).  Not  having  taken  into  account  these  losses, 
the  SNPSAM  model  generates  more  electric  power  and,  thus,  appears  more 
efficient  than  the  General  Electric  model. 

The  assumption  that  the  TEM  pump  does  not  act  like  a  heat  exchanger  (i.e., 
there  is  a  constant  temperature  boundary  condition  instead  of  a  heat  flow  boun¬ 
dary  condition)  is  only  a  rough  approximation.  General  Electric  reports  the  per¬ 
centage  of  heat  flow  from  the  primary  to  the  secondary  to  be  approximately  eight 
percent  of  the  rejected  heat.  Of  the  eight  percent  that  goes  to  the  secondary, 
seventeen  percent  is  rejected  by  the  TEM  pump  radiator,  or  1.4  percent  of  the 
total  rejected  power  (2).  The  1.4  percent  not  accounted  for  in  this  study  is 
justifiable  but  the  eight  percent  should  be  considered  in  future  studies. 
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The  A P  computed  for  the  secondary  heat  exchanger  was  done  so  with  the 
primary  heat  exchanger  hydraulic  model.  This  is  only  an  approximation  because 
the  fluid  properties  in  the  primary  and  the  secondary  heat  exchangers  will  be 
different  because  the  temperatures  are  different.  Since  the  friction  factor  depends 
on  the  properties  of  the  fluid,  diameter  of  piping,  and  type  of  material  used,  the 
pressure  loss  in  the  heat  exchangers  will  be  different.  It  is  recommended  that  a 
separate  heat  exchanger  hydraulic  model  be  added  to  the  revised  version  of 
SNPSAM. 

It  has  been  assumed  thus  far  that  General  Electric’s  results  are  a  suitable 
standard  of  comparison.  It  is  conceivable  that  some  of  the  disagreement  between 
results  of  SNPSAM  and  General  Electric  are  because  of  inaccuracies  in  the  General 
Electric  model. 

The  division  of  the  energy  conversion  assembly  at  the  plane  of  symmetry  is 
questionable  since  in  doing  so  the  fluid  that  flows  through  it  is  divided  in  half. 

The  plane  of  symmetry  halves  the  primary  heat  exchanger,  therefore,  its 
equivalent  diameter  is  altered.  In  this  study,  the  equivalent  diameter  of  the  pri¬ 
mary  heat  exchanger  was  computed  based  on  a  conduit  having  the  dimensions  of 
the  primary  heat  exchanger,  except  half  as  thick.  Since  the  heat  transported  from 
fluid  to  metal  is  a  function  of  the  heat  transfer  coefficient,  which  is  a  function  of 
the  equivalent  diameter,  the  choice  of  an  equivalent  diameter  must  be  justified. 
Further  investigation  into  the  choice  of  an  equivalent  diameter  is  warranted. 

It  is  worth  noting  that  the  unrevised  version  of  SNPSAM  bypasses  the  prob¬ 
lem  of  choosing  an  equivalent  diameter  by  not  dividing  the  system  about  its  mid¬ 
plane.  Instead  of  halving  the  system  and  doubling  the  results,  the  effective  surface 
area  of  the  heat  exchanger  was  doubled. 

The  secondary  coolant  model  developed  here  is  clearly  not  transient,  but  nei¬ 
ther  was  the  primary  coolant  model  described  in  Seo’s  thesis.  It  is  assumed  in 
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this  study  and  in  Seo’s  that  the  transient  behavior  outside  the  reactor  is  a  func¬ 
tion  of  the  thermoelectric  converters  only.  This  approach  was  selected  because  it 
allows  the  subroutines  to  be  written  independently  and  is  justified  because  heat  is 
transferred  more  quickly  from  the  fluid  than  from  the  thermoelectric  material. 

The  alternative  is  a  matrix  solution  of  the  entire  system.  Matrix  solution 
methods  are  used  in  the  reactor  model  and  in  the  thermoelectric  model.  These 
models  are  complicated  and  slow  in  generating  solutions.  It  is  estimated  that  the 
non-matrix  solution  method  (i.e.,  decoupling  the  system  models)  generates  results 
that  are  acceptable  for  non-rapid  transient  simulations  (3). 

To  better  insure  the  transient  behavior  of  the  fluid  does  not  come  into  play, 
Seo  let  the  minimum  time  step  allowed  in  a  simulation  be  equal  to  one  circulation 
time  of  the  coolant  fluid.  The  motive  for  doing  this  is  that  any  given  segment  of 
fluid  will  be  at  the  same  point  in  the  reactor,  and  encounter  the  same  conditions, 
as  it  did  one  circulation  time  in  the  past  (15).  This  approach  loses  its  meaning, 
however,  when  a  secondary  is  added  and  two  circulation  times  exist.  In  this 
study,  the  slowest  of  the  primary  and  secondary  circulation  times  is  taken  to  be 
the  minimum  time  step. 

A  deficiency  in  the  system  model  is  that  it  does  not  account  for  heat  losses 
that  are  caused  by  two  dimensional  heat  flow.  In  the  energy  conversion  assembly, 
for  example,  the  heat  lost  from  the  ends  and  edges  of  the  heat  exchanger  and  ther¬ 
mocouples  was  not  computed.  Also,  heat  loss  from  the  piping  is  not  considered. 
Both  SNPSAM  and  General  Electric  use  the  one  dimensional  approach. 

The  use  of  COMMON  statements  in  many  places  made  it  impossible  to  call 
certain  modules  (such  as  CONVEC,  HYDR,  DPIPE,  and  DPHEX;  section  3.1.2.3) 
for  use  in  the  secondary  module  without  altering  them.  Not  all,  but  many  of  the 
subroutines  would  be  made  more  accessible  by  making  better  use  of  arguments  in 
call  statements  (10). 
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Changes  that  were  made  to  the  input  file  and  to  the  source  code  often 
resulted  in  the  program  failing  because  the  thermocouple  model  could  not  con¬ 
verge  on  a  correct  solution  since  the  initial  conditions  were  out  of  range.  The 
method  of  making  changes  incrementally,  small  enough  so  that  SNPSAM  could 
converge  on  new  initial  conditions,  was  time-consuming.  It  should  be  possible  to 
write  a  command  file  that  executes  SNPSAM  for  cases  that  are  changed  incremen¬ 
tally,  eliminating  the  need  for  a  user  to  interactively  change  the  initial  conditions. 
Time  did  not  allow  this  to  be  done  in  this  study. 

The  transient  behavior  of  the  radiator  was  not  modeled  in  this  study.  It  is 
recommended  that  a  transient  model  for  the  radiator  be  developed.  Doing  so  will 
require  information  about  the  transient  behavior  of  heat  pipes  (such  information 
is  scarce).  It  will  also  require  that  a  two  dimensional  heat  transport  model  of  the 
radiator  be  written.  It  should  be  possible  to  generate  an  effective  thermoconduc¬ 
tivity  for  the  heat  pipes  and  treat  the  radiator  as  a  solid  with  that  thermoconduc¬ 
tivity. 

Although  SNPSAM  allows  flexibility  in  the  choice  of  dimensions  and  materi¬ 
als  used  in  the  space  reactor  design  being  modeled  (via  the  input  file),  it  does  not 
allow  for  changes  in  configurations  such  as  those  that  have  been  proposed  by  Gen¬ 
eral  Electric  since  Seo’s  work  and  since  this  study.  For  example,  since  this  study, 
General  Electric  has  chosen  to  stack  more  hot  and  cold  legs  together  in  the  energy 
conversion  assembly  (8).  It  is  therefore  recommended  that  SNPSAM  be  made 
more  modular  so  that  changes,  like  those  made  in  this  study,  can  be  effected  more 
easily. 

It  is  unlikely  that  SNPSAM  can  be  made  a  generic  computer  model  for  space 
reactors.  It  is  feasible,  however,  to  make  the  source  code  of  SNPSAM  more  modu¬ 
lar  and,  thus,  pliable  to  space  reactors  of  different  configurations.  Towards  this 
end  there  are  some  corrections  to  be  made  to  the  code  and  there  are  areas  that 
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require  further  study.  Based  on  the  conclusive  remarks  of  this  chapter,  a  sum¬ 
mary  of  recommendations  is  as  follows: 

■  Use  quadrature  basis  function  in  place  of  hat  basis  function  and  include 
higher  order  terms  in  finite  difference  scheme 

■  Write  a  subroutine  that  computes  the  secondary  head  loss  exclusively  (i.e.,  as 
opposed  to  using  the  subroutine  used  for  computing  the  head  loss  in  the  pri¬ 
mary) 

■  Write  in  heat  flow  boundary  condition  for  the  TEM  pump  (eliminating  the 
temperature  boundary  condition  used  in  this  study) 

■  Include  a  subroutine  that  accounts  for  the  heat  lost  between  channels  in  the 
radiator  (by  conduction  and  by  mass  transferred  through  thawing  holes) 

■  Determine  if  the  equivalent  diameter  used  in  the  heat  exchanger  is  appropri¬ 
ate 

■  Study  the  transient  behavior  of  the  radiator  and  include  a  subroutine  that 
accounts  for  this  behavior 

■  Implement  changes  that  make  more  use  of  arguments  in  call  statements  and 
find  other  ways  of  making  SNPSAM  more  modular 

■  Write  a  oommand  file  that  executes  SNPSAM  for  incrementally  varying 
cases,  aiding  in  the  finding  of  initial  conditions 

■  Add  documentation  to  the  code 
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Appendix  A:  SNPSAM  Input  File,  SYS.INP 


>>>>  SPACE  NUCLEAR  POWER  SYSTEM  SIMULATION  DATA  (Rev  GNE88M-.' 

6.8  MW,  Normal,  Transient, Li,  SiGe/GaP,  Heat  Flux  B.C., 
«**«  CALCULATION  OPTIONS  ###« 


Decay  Heat  option 
Operation  Option  ,- 

TE  Boundary’  Option 
Heat  Pipe  Sonic  Limit  ; 

Maanetic  Flux  Density  ,- 

MM* PLOTTING  INTERVAL  *### 
Printing  intervals 
#*M  FLUID  AND  SOLID  M*i 
Coolant  ID  Number 
Pump  Duct  Material  ; 

Heat  Exchanger  Wall  Mat.  ,- 

Compliant  Pad  Material  ,- 

MM  CORE  PARAMETERS  FOR  KINET: 
»»  REACTIVITY  *» 


0  1  — >  decay  heat,  other 

0  0  —  >  steady  state,  1  — 

0  0  -->  heat  flux,  1  — 

0  0  — >  off,  1  — 

0  0  — >  Temp.  Depend.  1  — 

1 

1 

2 

2 

1 

CS  AND  THERMAL  CALCULATION  M* 


Type  of  Reactivity  Function;  2 
Number  of  changes  in  react.  ,-  2 
Time  Reactivity 

1.0  0.0 

2.0  -1.0e-4 


Simulation  End  Time (sec)  ,- 
**  INITIAL  CONDITIONS  »» 
Initial  Power (W)  ,- 

Core  Inlet  Temperature ( K ) 
Fuel  Average  Temperature (K) ; 
Cladding  Average  Temp.(K)  ,- 
Mass  Flow  Rate (kg/sec)  ,- 

secondary  Inlet  Temp.(K) 
Secondary  Outlet  Temp . (K) 
Secondary  Mass  Flow( kg/sec); 


250.0 

6 . 8000e6 
1217.775 
1654.724 
1353.489 
12.08629 
736.9515 
805.2298 
21.92900 


**  CORE  DESIGN  PARAMETERS  ** 

Total  Heat  Transfer  F->CL  ;  2.167e4 
"  CL->C0  ;  9. 531e4 

Fuel  Specific  Heat(J/kgK)  ;1.465e2 
Cladding  ”  ; 0 . 276 

Coolant  volume  in  Core ( m3 )  ;0.015 
Total  Cladding  Weight  (kg)  ;2.55el 
Total  Fuel  Weight  (kc)  ,-1.09d2 

•’  Reactivity  Temperature  Coefficient  ** 

Doppler  Reactivity  Coeff.  ,-2.4  7e-7 
Core  Expansion  Coeff.  ;-1.26e-5 

Coolant  Expansion  Coeff.  3.25e-6 
**  Time  Step  for  Integration  *» 

Maximum  Timestep  Allowed  ;  5.0e0 
Minimum  "  "  ;  1.0e-7 

Delta  t  ;  0.01 

1-Fixed,  2-Variable  Time  Sd;  2 
##M  CORE  PARAMETERS  FOR  PRESSURE  DROP  CALCULATION  #### 


Channel  Height (m) 

Wire  Wrap  Lead(m) 

Fuel  Rod  pitch(m) 

Fuel  Rod  Diameter(m) 
Control  Rod  FTF  Length (m) 
Rx. Vessel  Inner  Dia.(m) 
Equiv.  Baffle  Dia.(m) 
Number  of  Fuel  Rods 
#1##  EM  PUMP  DATA  FOR  PUMP 
Duct  Height (m) 


;0 . 454 
;0. 37 
,•0.007981 
;0. 0074 
,-0.0465 
;0.358 
;0.328 
,-  1296 

HEAD  CALCULATION  M## 
,-0.0268 


A.l 


)  <<<< 


->  normal 
transient 
temperature 
on 

constant 


Duct  Width (m) 

Duct  Length(m) 

Duct  Wall  Thickness (m) 
Magnetic  Flux  at  To (gauss) 
Temperature  To 
Number  of  TE  couples 


0.0178 

0.2550 

0.0006 

1000.0 

950.0 

50 


f *M  EM  PUMP  DATA  FCR  SECONDARY  PUMP  HEAD  CALCULATION  4444 
Duct  Height (m>  ,-0.0171 

Duct  Width (m)  ; 0.0116 

Duct  Length(m)  ,-0.2550 

Duct  Wall  Thickness (m)  ,-0.00064 

Magnetic  Flux  at  To(gauss)  ,-1000.0 
Temperature  To  ,-  950.0 

Number  of  TE  couples  ,-  50 

4444  PERMANENT  MAGNET  DATA  FOR  TEMPERATURE  DISTRIBUTION  4### 


Number  of  Multifoils 
Emissivity  of  multi foils 
Foil  Thickness 
Gap  width  between  foils 
1/2  Width  of  Magnet 
1/2  Height  of  Magnet 
Thermal  Cond.  of  Magnet 
Number  of  meshes  in  Magnet 
Number  of  meshes  in  Foil 
Number  of  meshes  i  height 
Number  of  Terms  in  Sol. 


30.0 

0.4 

0.001 

0.004 

0.5 

5.80 

0.4 

5 

2 

5 

25 


*44*  HEAT  EXCHANGER  DATA  FOR  THERMAL  CALCULATION  4*44 
Channel  Inner  Width(m)  ,0.2133 

Channel  Outer  Width(m)  ,-0.2138 

Channel  Inner  Height(m)  ,-0.004212 
Channel  Length(m)  ,-2.7736 

Number  of  channel (f 10. 5)  ;  1.0 

Wall  Thickness (m)  ,-0.0005 

4  of  segment/channel  ,-10 

*  of  TE  couple/segment  ,-  2662 

4444  HEAT  EXCHANGER  DATA  FOR  PRESSURE  DROP  CALCULATION  444# 
Header  length(m)  ,-0.2133 

Header  inlet  diameter(m)  ,-0.0413 

Header  end  diameter(m)  ,0.019 

Channel  Length(m)  ,-2.7736 

4444  PIPS  DATA  FCR  PRESSURE  DROP  CALCULATION  ### 4 

41  Pipe  length(m)  ,-1.2 

42  Pipe  length(m)  ,-0.3 

43  Pipe  length(m)  ,-0.3 

44  Pipe  length(m)  ,1.2 

41  Pipe  diameter (m)  ,-0.0826 

*2  Pipe  diameter (m)  ,-0.0413 

43  Pipe  diameter(m)  ,-0.0413 

44  Pipe  diameter(m)  ,-0.0826 

44**4  RADIATOR  DATA  FOR  THERMAL  CALCULATION  44*4 
Channel  Inner  Width<m)  ,0.1500 

Channel  Outer  Wldth(m)  ,0.1505 

Channel  Inner  Height(m)  ,0.0090 

Channel  Length(m)  ,-25.0 

Number  of  channel (f  10. 5)  ,-  1.0 

Wall  Thickness (m)  ,-0.0005 

4  of  segment/channel  ,-10 

4  of  TE  couple/segment  ,-  2662 

4*44  SECONDARY  HEAT  EXCHANGER  DATA  FCR  PRESSURE  DROP  CALCULA 
Header  length(m)  ,0.8 


Header  inlet  diameterfm)  ,-0.015 

Header  end  diameter (m)  , -0.015 

Channel  Length(m)  ;10.0 

«M4  PIPE  DATA  FOR  SECONDARY  PRESSURE  DROP  CALCULATION 
*1  Pipe  length<m)  ;1.0 

*2  Pipe  length(m)  ,25.0 

♦  3  Pipe  length(m)  ,-25.0 

»4  Pipe  length(m)  ;1.0 

♦1  Pipe  diameter(m)  ,-0.030 

#2  Pipe  diameter(m)  ,-0.030 

#3  Pipe  diameter(m)  ,-0.030 

*4  Pipe  diameter(m)  ,-0.030 

*♦4#  TE  DATA  FOR  ENERGY  CONVERSION  CALCULATION  ##*# 

***  SiGe/GaP  Semiconductor 

Option  is  i  2  — >  1  -  fixed  Tc,  2  -  radiative 

TE  material  ;  1  — >  1  -  SiGe/GaP,  2  -  SiGe 

Number  of  Element  ;10 

"  Node  ,-11 

Boundary  Condition  ***** 

Left  Boundary  Flag  ;  0 

•Do  you  have  external  heat?  ;  1  -**>  0  -  no,  1  -  yes 

Number  of  changes  ,-  1 

put  end  time  t(i),  coefficient  a(i)  and  b(i)  in  Ax  +  B. 

2000.0  0.0  0.0 
*****  Length  of  Element  ***** 

0.038  0.038  0.038  0.038  0.038  0.038  0.038  0.038  0.038  0.038 

**»»  x  coordinat  of  nodes  **** 

0.00  0.038  0.076  0.114  0.152  0.190  0.228  0.266  0.304  0.342 

0.38 

****  Initial  Condition  **** 

N-LEG; 

1311.0851266.8851221.9031176.1391129.5911082.2601034.145  985.245  935.559  885.088 
833.830 
P-LEG; 

1311.5631267.0431221.8011175.8381129.1511081.7411033.608  984.750  935.168  884.862 
833.830 
N-LEG, 

1298.2111254.5071210.0321164.7851118.7674071.9751024.410  976.071  926.958  877.070 
826.406 
P-LEG , 

1298.6451254.6021209.8531164.3981118.2361071.3681023.792  975. 50e  926.516  876.815 
826.406 
N-LEG; 

1284.8981241.5771197.4901152.6361107.0141060.6251013.468  965.542  916.847  867.381 
817.146 
P-LEG; 

1285.3091241.6221197.2411152.1651106.3941059.9271012.764  964.905  916.349  86"\096 
817.146 
N-LEG; 

1271.3091228.3861184.7031140.2601095.0561049.0911002.364  954.875  906.623  857.607 
807.827 
P-LEG; 

1271.6941228.3811184.3861139.7091094.3481048.3051001.578  954.167  906.071  857.292 
807.827 
N-LEG ; 

1258 . 8891216 . 3091172^  9731128.8831084.0361038.434  992.075  944.959  897.085  848.453 
799.062 
P-LEG; 

1259.2511216.2561172.5911128.2541083.2451037.564  991.210  944.183  896.483  848.109 
799.062 


N-LEG; 

1246.6371204 . 3681161.3501117.5841073 . 0681027 . 803 
790.215 
P-LEG ; 

1246.0221203.4121160.1431116.2151071.6261026.377 

'90.215 

N-LEG; 

1233 . 6501191.7601149. 1291105 . 7591061.6471016 . 794 
781.391 
P-LEG ; 

1233.0771190.8211147.9181104.3691060.1721015.329 

781.391 

N-LEG; 

1221. 9671180 .4241138. 1491095. 1411051. 4001006 . 926 
773.537 


981.788 
980. 4e7 

971.199 

969.837 

961.713 


935.022 

933.896 

924.861 

923.698 

915.775 


887.505 

886 . 664 

877.781 

876.910 

869.098 


839.236 

838. *70 

829.958 

829.475 

921.685 


P-LEG; 


1221.4361179.5041136.9381093.7371049.9011005.430 

773.537 

N-LEG; 

1209.0781168.0311126.2461083.7261040.475  996.492 


765 . 665 


960.323  914.581 


951.780  906.341 


P-LEG ; 

1208 . 6131167 . 1471125 . 0481082.3171038 . 958 
765.665 
N-LEG; 

1197.2231157.0541116.1971074.6521032.419 

764.533 

P-LEG; 

1196.8791156.2921115.1181073.3551031.004 

764.533 

**”  cross  sectional  area  ******* 


Area  of  the  lea.  An : Ap :  ;  0.104 

*******  Radiator  parameters  ******** 
Surface  Emissivity  ;  0.85 

Ambient  Temperature  ;  250.0 

Shape  Factor  ,-  0.95 

Area  Ratio  ;  30.82 

*******  External  Load  ***»»•** 


994.971  950.358  905.119 
989.496  945.885  901.583 
988.065  944.537  900.420 

0.104 


868.202  821.188 

860.174  813.282 

859.257  812.772 

856.591  810.907 

855.714  810.419 


Number  of  changes  ,-  1 

put  end  time  t(i),  coefficient  a(i)  and  b(i)  in  Ax 


1000.0  0.0  0.0161 

14*4  TE  DATA  FOR  EM  POMP  CALCULATION  4444 
**  SiGe/GaP  DATA  ’* 


Option  is 
TE  material 
Number  of  Element 
"  Node 


;  2  —  >  1 

;  1  — >  1 

i  B 

;  9 


*****  Boundary  Condition  ***** 

Left  Boundary  Flag  ;  1 

•Do  you  have  external  heat?  ;  1  — >  0  - 

Number  of  changes  ;  1 

put  end  time  t(i),  coefficient  a(i)  and 

2000.0  0.0  0.0 


fixed  Tc, 
SiGe/GaP, 


no,  1  -  yes 
b(i)  in  Ax 


B. 


2  - 
2  - 


B. 


radiative 

SiGe 


Length  of  Element  ***** 


0.080  0.080 
****  x  coordinat 

0.080 
of  nodes 

0.080 

*  *  *  » 

0.080 

0.080 

0.080 

0.080 

0.00  0.080 

D.160 

0.240 

0.320 

0.400 

0.480 

0.560 

****  Initial  Condition  »»** 

N-LEG ; 

1352. 2201308. 414 1255. 3421194. 0401125. 0851C48. 672  964.64'  872.493 


0.640 


A.  I 


8 


% 


?-LEG ; 

.352. 2201306 . 8881252. 9451191. 2041122. 121104  5. 852  962. 276  870.989 
******  Radiator  parameters  »*•«»»»•« 

Surface  Emissivity  ,-  0.85 

Ambient  Temperature  ,-  250.0 

Shape  Factor  ,-  0.95 

Area  ratio  ,- 15 . 0 

.......  pressure  drop  weighting  factor  .......... 

Weighting  Factor  for  HYDR  ,-0.9 

Weighting  Factor  for  SCNDRY  ;0.9 

.......  Maximum  Temperatures  ........... 

Fuel  Melting  Point  / 3 073.0 

Cladding  Melting  Point  ; 2370.0 

Coolant  Boiling  Point  ,-2590.0 

TE(SiGe/GaP)  Melting  Point  ,-2500.0 
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Appendix  B:  SNPSAM  Output  File,  PLOT.TMP* 


0.000 
0 . lOOOOOOD-35 
0. 1350366D-04 
0 . 6261468D-04 
0.7971588D+02 
0.12173850+02 
0.6326531D+01 
0.7373537D+03 
0.000 
0. 1000000D-35 
0 . 13503660+04 
0.6261468D+04 
0 . 7971588D+02 
0 . 12173850+02 
0. 63265310*01 
0 . 7373537D+03 
3.400 
0.7000000D-05 
0 . 13493050*04 
0.6247194D+04 
0.82713950+02 
0. 12147670+02 
0.6315481D+01 
0 . 7369204D+03 
6.805 
0.2402500D-04 
0 . 1349996D+04 
0.62425510-04 
0 . 8290002D+02 
0.12152830-02 
0 . 6322364D-01 
0.73727170+03 
10.210 
0.4105000D-04 
0 . 1350322D+04 
0 . 6247305D+04 
0. 8313278D+02 
0. 1215381D+02 
0. 6325237D+01 
0 . 7373335D+03 
13.620 
0.5810000D-04 
0. 1350835D+04 
0. 6249533D+04 
0 . 8328609D+02 
0 . 1215992D+02 
0.6331706D+01 
0 . 7374875D+03 
17.020 
0. 7510000D-04 
0.13513620+04 
0.62553890+04 
0 . 8325962D+02 
0. 12166650+02 
0 . 6330848D+01 
0.73761990+03 


0 . 1000000D-3  5 

0 . 1285010D-04 

0 . 6441014D-Q4 

0 . 3774698D-02 

0 . 6714463D-04 
0. 1018572D-00 

0 . 80433580-03 

O.OOOOOOOD-OO 
0. 1654748D-04 
0.13522450-04 

0 . 1081213D+05 

0 . 4148327D-03 

0 . 1246422D-04 

0 . 3725466D+01 

0 . 6800Q00D-04 

0 . 1353502D-04 
0. 77115460-03 
0. 1729089D-05 
0.6060878D-04 
0. 7897633D-03 
0.7630722D-03 

0.1217575D+04 
0. 1795466D-03 
0 . “061~18D-02 
0. 3287700D-05 
0. 6171401D+01 
0 . 1171445D-01 

0 . 1000000D-35 
0. 1285010D+04 
0. 6441014D-04 

0 . 3774698D+02 
0. 6714463D-04 

0 . 1018572D-00 
0.80433580+03 

0. OOOOOOOD+OO 

0 . 1654748D-04 
0. 1352245D+04 

0 . 1081213D+05 

0 . 4148327D+03 

0 . 1246422D+04 
0.3725466D+01 

0. 6800000D-04 

0 . 1353502D-04 

0 . 7711546D-03 
0. 17290890-05 
0. 60608780-04 
0.7897633D-03 

0 . 7630722D-03 

0.12175750+04 
0.17954660-03 
0 . 70617180-02 
0 . 3287700D-05 
0.61714010+01 
0.11714450-01 

0.9778266D-05 

0.12844160+04 

0 . 6426337D+04 
0.37661510-02 
0.6697123D+04 
0. 1016792D+00 

0 . 8050664D+03 

0.2778266D-05 
0.16546310-04 
0.13512570-04 
0.1078720D-05 
0. 4145678D-03 

0 . 1245936D+04 
0.3659411D-01 

0.6812353D-04 

0.13530100-04 

0.77091090-03 

0 . 1725123D-05 
0.58989880-04 

0. 7897633D-03 
0.76322'’  .0*73 

0 . 121~826D-04 
0. 1* 914290-03 
0. 7028845D-02 
0.32750270+05 
0. 6183571D+01 
0.1115853D-01 

0.13488280-04 

0 . 1284890D-04 
0. 6421551D-04 

0 . 3763220D-02 
0.6709195D-04 
0.10179010-00 
0.80514950-03 

-0.10536720-04 

0 . 1655615D-04 

0 . 1351953D-04 

0 . 1077907D+05 

0 . 4152647D+03 
0.1246435D-04 
0.3677341D+01 

0 . 6821053D-04 

0. 1353571D+04 
0.7710677D-03 

0 . 1723825D-05 

0 . 5710420D-04 

0. 7897633D-03 
0.7634447D-03 

0. 1217992D-04 
C. 1790002D-03 

0 . 7036198D-02 
0.32777730-05 
0.61827020+01 
0.1116697D-01 

0 . 1854948D-04 
0.12851380+04 
0.6426447D+04 
0. 3766174D+02 

0 . 6712742D+04 
0. 1018363D+00 

0 . 8052783D+03 

-0.22500520-04 

0. 1656501D-04 

0 . 1352284D-04 
0.1078754D+05 
0. 415694 0D-03 

0 . 1246631D+04 
0.3675525D+01 

0.6832752D-04 

0 . 1353929D-04 
0.77127800-03 

0 . 1725165D+05 

0. 5573868D+04 
0.7897633D-03 
0.7635344D-03 

0 . 1218305D-04 

0 . 1791417D-03 
0. 7037079D-02 
0 . 3278286D+05 
0.61858320+01 

0 . 1113953D-01 

0.2156312D-04 
0. 1285553D+04 
0. 6428740D+04 
0.3767551D+02 

0 . 6719173D+04 
0. 1019405D+00 
0. 8053922D+03 

-0 . 3653688D-04 

0 . 1657508D+04 

0. 1352800D+04 
0.1079155D+05 

0 . 4163578D+03 

0 . 1247076D+04 

0. 36801620+01 

0.6843407D+04 

0. 1354458D+04 
0.7712780D+03 
0.17257960-05 
0.58875560-04 

0 . 7898793D-03 

0. 7636648D-03 

0.12187020+04 

0.1792071D+03 

0.7043814D+02 

0 . 3281389D+05 
0. 6189721D+01 
0. 1114339D-01 

0.23911450-04 

0. 1286014D+04 

0 . 6434769D+04 

0. 3771174D+02 
0.6724280D+04 

0 . 1019267D+00 
0.8055262D-03 

-0 . 5118855D-04 

0. 1658552D+04 

0 . 1353326D+04 

0 . 1080195D+05 
0.41689260+03 
0.12471140*04 

0 . 3682638D+01 

0.68538320-04 

0. 1355021D-04 

0 . 7713676D-03 

0. 1727443D-05 

0 . 6069297D-04 

0. 7900028D-03 
0.76379230*03 

0.1219195D+04 

0.1793805D+03 

0 . 7050880D+02 

0 . 3284763D-05 
0.61928960-01 
0.11171530-01 

•  Output  parameters  for  PLOT.TMP  are  listed  in  Appendix  F. 


Appendix  C:  Subroutine  SECONDARY 


3 

14 

15 


39 

:o 


13 

34 


cccocccccccccccccccccccccccccccc 
c  SUBROUTINE  SECONDARY 


'.u,  reyf ,  prar.:  . 


ccccccccccccccccccccccccccccccccccccccccccc 
subroutine  secondary (bcrnn) 
implicit  real*8  (a-h,o-:) 

real '8  mdco,mvf ,mul,muv 
dimension  otfs( 30 ) ,qdct ; 30 ) 
common/f sid/if luid. idl , isol , iso2 
common/hxam/tid , tod, thi, tlen. tn, tv, ntep 
»  pechtf , shia . Dit . dea, anu , nsa 

COMMON/NMBR/MENU ,  NE ,  NN,  NS  .  ISTEP ,  STEP ,  PSTEP  ,  "CUNT ,  ITEM 
common/tmpr/ttn( 30,21), tt?( 30 , 21) , ttf ( 3 , 30) , ttv( 3,30), 

*  tcrr;30) 

COMMON  / C0RTP5/  TECA , TECA0 , TECI , TECX . TEFA , TEFA0 , TELA . TE1 
CMF , 2ML , VC , ULA , UCA , SOCA , Zc 1C0R , ZRCCX , ZPCCRA , ZEGXC , ZTEC3L . 
CX. ZFLCR0 

common/ scndry/TES I , TESX, iscnary , zrlsec , hfsec . ttns(  30) 

□ran ; cp , vis , tk ) -cp’vis/tk 

rev ( vel , dea , den , vis ) -deq’ vel ’den/vis 

pecht ( deq , vel , den , cp , tk ) -deq’vei *den ’cp / tk 

mdct-ZFLSEC/12 . 0/tn/2 . 


ZVC 


26 

999 

tcheck-TESI 

call  radiator 

28 

tin-TESI 

29 

c 

30 

ffa«tid*thi 

31 

hper-2 . 0*tid 

3  2 

hta-1 . 0’tlen/r.sg’tid 

33 

c 

34 

tmpdr-1 . 0-  ( TES I -TESX )  /ti.n/nsg 

35 

tout-tin’tmpdr 

36 

taf-( tin* tout) /2. 0 

3" 

c 

38 

do  300  i-l,nsg 

39 

icc-1 

40 

c 

41 

call  fluid ( i fluid, taf , pv.mwf , rhoi 

42 

*tkf ,cpf ,resf ) 

43 

c 

44 

qdothf -hfsec *ntep 

45 

c 

46 

if(i.eq.  {nsg-l-iscndry )  )the.n 

4* 

deq-4 . 0’f fa/hcer 

48 

veif-mdot/( f fa’rholf) 

49 

reyf-rey ( vel f , deq , rhol f 

80 

pranf-pran(cp: ,mul,tkf ) 

51 

pecht  f  -  pecht  ( dec ,  vel  f ,  r; 

52 

call  CONVEC 

53 

hf-anu’tkf/deq 

54 

tw-taf*qdothf ,'( hr ’hta ) 

55 

c 

56 

tva-tv 

5' 

:60 

call  soiidf isol, tva. amp 

58 

twl-tw- 'qdothf * ' tod-tid 

59 

twar.-(tv-twl)  2.3 

60 

i f  ( dabs ;  tvan-  tva  ; ,  tva .  g1 

nov , 


rhos ,  cps  ,  t.<s  ,  tec ,  rec 


C.l 


Con  avafiatte  to  DTIC  does  not 
permit  full?  legible  reproduction 


r  n 


61 

62 

63 

64 
6  3 
66 

67 

68 
69 
'0 


'3 

'4 


~6 


78 
-  9 
8C 
81 
82 
S3 

84 

85 

86 
87 
98 

89 

90 

91 

92 

93 

94 

95 


c 

261 


c 


£va-twan 
go  to  260 
end  if 

twa-twl 

call  solid( isol,  tva . a.tp , rhos , cps ,  tics .  tec,  res ) 
ttns  ( i ; -twl-  ( qdothf  *  0 . 2  0 1 2  5 )  /  ( tks  *  hta ) 
twan-(twl-ttns(i) )/2. 0 
if  (dabs; tw&n-tvai/tva . gt . 3 . 301) then 
tva- tva  n 
go  to  261 
end  if 

bcrnn-ttns(i) 


end  if 
c 

do  866  j-nsg~l-iscndry , 1 , -1 
966  qdot(j ) -qdothf 
c 

tout-tin-qdct ( i ) /cp  f /mdot 
taf- ( tout-tin ) /2 . 
tin-tout 
ttfs(i)-tin 
c 

800  continue 
c 

TESX-tout 

reierror-acs ( {tcheck-TES:;,  TES2; 
c  write (•»,*)  relerror 

if(reierror.gt.3,005)  go  to  999 
call  hydrsec 
return 
end 
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Appendix  D:  Subroutine  RADIATOR 


3 

4 

5 

6 

8 

9 

10 


13 

14 

15 
1.6 

18 

19 

CO 

21 

22 

23 

24 

25 

26 
27 
23 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 
5“ 

58 

59 

60 


cccccccccoccccccccoccccccocccccrcoccccccor:: 
c  c 

c  SUBROUTINE  RADIATOR  c 

c  c 

cccccccccccccccccccccccccccccccccccccccccco 
subroutine  radiator 
implicit  real’8  (a-h,o-z) 
real *8  mdct.mwf ,mul,muv 
dimension  ttnr;20) 
common/rdtr  emss , tsur , sha? 
ccmmon/zero/heatfx, const 

conunon/fsid/ifluid,  idl ,  iscl ,  iso2 
common/hxdmrad/tidrad .  todrad, thirac .  tlenrad,  tr.raa.  tvrad ,  r. 
*  inu , rey f . tran  f , pecht r , shid , pit , dec ,  anu , nsa 

COMMON/NMBR/MENU , ME , NN . NS , ISTE? , STEP , PSTE? . XC0UNT , ITEM- 
CO-TON  /COP. TP 5/  TECA ,  TECA0  .  TEC!  .  TECX .  TEFA ,  TEE A0  .  TELA .  IS 
CMF . CML. VC , ULA , UCA , ROCA. OFLCOR , IRC CM . IPCCRA . IE3XC , ITECBL 
CX.IFLCR0 

common /scndry/TES I, TESM, iscndry, zrlsec, hfsec, tcns( 30 ; 
common/radiator/tradav 
external  cold 
c 

pran ( c? , vis , tk ) -cp’vis/tk 
rey ( vel ,  dec ,  den ,  vis )  -aeq’vei’den/vis 
pecht ( deq ,  vel , den , cp , tk ) - dec* vel ’den • cp/tk 
c 

tn-tnrad 
nsg-nsgrad 
tid-tidrad 
tod-todrad 
thi-thirad 
tlen-tlenrad 
ntep-nteprad 
c 

mdot-ZFLSEC/12 . 0/tn/2 . 
tin-TESX 

tmpdr-1 . 0-(TESX-TESI)/tir./nsg 
c 

ffa-tid*thi 
hper-2.0*tid 

c  hta-2 . 0*tlen/nsg*tid 

hta-1 . 0’tlen/ nsg’tid 
c 

do  800  i-l,nsg 
icc-1 
c 

tout-tin’tmpdr 
taf-( tin-tout )/2.0 


teprad, 

rad 

0.CF.CL.Q 
VOLFC , IHC 


100 

c 


call  fluid ( if luid,  taf ,  pv.mwf , rholf  ,muv, mul , hr'c,  sio,rk, rhov, 
♦tkf.cpf ,resf) 

qdot-mdot*cpf ’ ( tin-tout ) 
qdotl-qdot/ntep 

Heatfx-CDCT/ntep 

tcold-zercin( tsur, 3000 . 0 , coid ,  I . Oe-9 ) 
write;*.’)  tcold 
ttnr ( i )-tcold 


I 


i 


D.l 


61 

tradav-(ttnr;l)*ttnr;  r.sg;  ;/2. 

62 

C 

* 

63 

tva-ttnr(i) 

64 

200 

call  solid(  isol .  :va  .  amp .  rncs ,  ops ,  tics  . 

63 

tv--tnr  ( i )  -  ( qdoc  *  0.00135)  ( tks  *  hta ) 

66 

tvan-  ( tv-ttnr  ( i ) )  /2 . 0 

67 

if ( dabs (tvan-tva) /tva. gt. 0.001)  then 

66 

tva-tvan 

69 

go  to  200 

70 

endif 

7  X 

C 

'2 

tva- tv 

"3 

tVl-tW 

74 

250 

call  soiid(  iso2 ,  tva ,  amp .  rhos .  cps .  tics . 

*5 

tv-twl-  ( qdot’  ( tcd-tid )';  /( txs*hta ) 

'6 

tvan- ( tv- twl ) /2 . 0 

7  7 

if ; dabs (twan-tva; /tva. gt. 0. 001)  then 

"8 

tva -tvan 

79 

go  to  250 

80 

er.d  if 

81 

C 

82 

dec- 4 .  0*ffa/hper 

83 

velf-mdot/( ffa-rholf ) 

84 

c 

85 

rev f- rev ( velf , deq , rholf . mul ; 

86 

pra.nf-prar.  ( cpf ,  mul ,  tkf ; 

87 

pechtf-pecht ( deq , velf , rholf .cpf, tkf; 

88 

c 

89 

call  CCXVECRAD 

90 

c 

91 

hf-anu-tkf/deq 

92 

c 

93 

tcutn-( (maot’cpf  -  hr *hta, 2. 0 ) -tin  - 

94 

*  Cmdot'cpf  -  hf-hta/2.0) 

95 

c 

96 

c 

vrite  ( * ,  - )  tout ,  tout.-. 

97 

aelt-toutn-tout 

98 

i f ;  dabs ( toutn- tout ) /tout . gt . 0 . 000  5 )  t 

99 

if(icc.gt.lCO)  then 

100 

c 

if (icc.eq.  101)  print  *,  '  RADIAT 

101 

abstt-dabs; toutn-tcut ;/tout 

102 

c 

print  ♦,  '  Difference  «',  abstt, 

103 

tout-tout-0 . 0057*delt 

104 

if ( icc.gt. 300 )  then 

105 

print  *,  'Failed  Again  !!!!' 

106 

if(abstt. It. 0.001;  go  to  ”9~ 

107 

stop 

108 

end  if 

109 

go  to  798 

110 

endif 

111 

c 

112 

if (dabs (toutn- tout),  tout. It. 0 . 002 ; 

113 

tout-tout-0 . C053*delt 

114 

else 

115 

tout-tout-0 . 0093’delt 

116 

endif 

117 

c 

118 

798 

icc-icc-1 

119 

go  to  100 

120 

endif 

TvoSIobl.  w  0«c  a"’“l 


121  c 

122  trodr-tout/tin 

123  3 

124  tin-tout 

125  c  write; »,*;  i.tin 

126  c 

127  800  continue 

128  c 

129  TESI-tOUt 

130  c 

131  end 
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Appendix  E:  Subroutine  EMPMP 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 

25 

26 
“4  — 

■ 

29 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c  c 

c  SUBROUTINE  EMPMP  (Rev  GNF88M-X)  c 

c  c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC'CCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCC:: 

c 

subroutine  EMPMP 
implicit  real*8  (a-h,o-z) 
real *8  mdot 

real* 8  mdotsec 


ccmmon/option/il,12,i3,i4 . i5,i6,i7,i8,i9,il0, 

*  zl,z2.z3,z4,z5,z6,z7,z8,z9,z!0 
common/f sid/if luid, isolid, id2,id3 
common/empmp/a , b , tl , wt , bf i eld , nte , phead , ef  f a 

comnvcn/empmpsec/asec,bsec  tisec ,vtsec,bfieldsec, ntesec. pheadsec , ef 
common/magl/bstand,  tstand.  foil,  emm,  t)an,  tmag,  trad,  to 
ccmmon/mag2/amg,bmg,  cmg,  n>  1 ,  nx2,nymg,nsmg,  foilth,gapth 
common/ti t/temp 

common/hydr/rv, rf , rbar , crro, rhol , weight 

CCMMON/GLOBP/GAMN (21,21;  ,GAMP(  21, 21) ,  GCMN(  21 , 21 )  ,GCMP '  21 , 21 ) , 

*  GBVN ( 21 ) , GBi? ( 21 ) ,  TVN( 3 , 21) , 7VP( 3 , 21 ) , 2CM( 21 , 3 ) , 

*  GAMTN( 21 , 21 ; , GAMTP (21,21), GBVTN ( 21 ) , G3VTP ( 21 ) , 

*  GAMG (21,21),  G3VG  (  21 )  ,  S0LV{  21 ) 

COMMCN/NMBRP/MENU, NE,NN, NS, ISTEP, STEP, PSTEP, IC0UNT, ITEM 
COMMON/BNDYP/BC0L<  3  )  ,  BCC  R(  3  ) ,  3CTL(  3  ) ,  BCTR  (  3  ) ,  I8CL ,  I3CR 
CCMMGN/PHYSP/AREAN , AREAf , THCCN , THC0P , RESTN , RESTP , CAPAN , CAPA? , 

*  ALPHA( 3,90), NCDT , RST (3,90), CON (3,90) 
COMMON/RDTRP/EMSS , TSUR , i HAP , FACTOR 

CCMMON/GNRLP/REXT , CRR , TAUN , TAUP , BCRN , BCRP , XV ( 2 1 ) , HV ( 2 I ) 
CCMMON/RSLTP/EFFV ( 3 ) , PWEV( 3 ) , CRRV( 3 ) , V0LV( 3 ) , PWIY( 3 ) 

*  , PWRV ( 3 ) 

CCMM0N/2ER0P/HF , CONST 

COMMON  /C0RTP5/  TECA , TECAC , TEC! , TECX , TEFA , TEFA0 , TELA . TELAO , CF , CL , Q 
CMF ,  QML ,  VC ,  ULA ,  UCA ,  ROCA ,  ZFI.CCR ,  ZRCCX ,  ZPCORA ,  ZEGXC ,  ZTEC3L ,  1YOLFC ,  ZHC 
CX, ZFLCRO 

ccmmon/scndry/TESI , TESX , ircadry , zf lsec , hf sec , ttns ( 30 ) 
save  cpf,t00,re 
c 

mdot-ZFLCOR/12. /2. 

mdotsec-ZFLSEC/12 . /2 . 
tO-TECX-t-273.0 
tOsec-TESX 
c 

if(ilO.eq.l)  then 
tOO-TECX 

call  solid( isolid , to , dl , d2 , d3 , d4 , d5 , ress ) 

call  f luid( i fluid , to , dl , d2 , rhol , d3 , d4 , d5 , d6 , d7 , d8 , d9 , cpf , res f ) 

call  solid( isolid, tOsec, dl,d2 ,d3,d4 ,d5 , resssec) 

call  fluid< if luid, tOsec, dl,d2,rholsee,d3,d4,d5,d6,d7,d8,d9, cpf ,r 

*  resfsec) 
c 

re-(resf*a/(b*tl) )/2.  +  resfsec*asec/(bsec*tlsec; 
rw-(ress*a/(2.*tl»vrt) )/: .  *  resssec*asec/( 2 . *tlsec*wtsec) 
rf-(10,0»ress*a/(2. *tl*vt) )/2.  «■  10.0'resssec’asec/ 

*  (2. *tlsec*wtsec) 
rbar- rw* re- r f * re-rv* r f 
ref f-re*rv*rf /rbar 
reff-reff*2.0 

c  re-resf*a/(b*tl) 
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atoflatto  to  OTIC 


c  r*»-ress*a/(2.  *tl*wt) 

c  rf-10 . 0*re'ss»a/{  2. ’tl’v-t; 

c  rbar-rw*re»ri»re*TW’rf 

c  ref f-re*rv*rf/rbar 

c  reff-reff *2. 0 

c 

call  TEP(tO,reff) 
c 

crro-crrv( istep ) «nte 
trad- f/n ( istep- 1 ,  nn ) 
c 

if;i8.eq.l)  then 
bfield-38 
go  to  500 
endif 
c 

call  magtem 
c 

500  continue 
endif 
c 

q-mdot/rhol 

qsec-mdotsec/rholsec 

phead-bfield/(  (b/2.-bsec;,rbar',1.0e4)  *  (rw»rfTcrro  - 

*  bfield/(1.0e4*(b/:..-bsec)  )*(rw-rf)’,(q/2.-asec) ) 

ef  f-phead*  { q/2 .  -qsec  5/crrC'/' ' fcf ield*  ( a/2 .  -qsec ; /I .  e4/(  b/'2 .  -bsec ;  - 

*  * l.e4»(b/2.-bsec)*re*phead/bfield) 

c  pnead-bfield/(b’rbar’1.0e-4)  *  (rv*rf*crro  - 

c  -  bfield,'(l.  Qe4*b)  •  (rv-rf )  *q) 

c  eff-phead*a/crro/(bfield*q/l.e4/b  -  l.e4*b* 

c  ’  re'pnead/bf ield ) 

ef  f  a-ef  f  « ef fv ( istep ) 
c 

phead-phead*2.0 

c 

en>pdt-pwiv(  istep)  *nte/(mdc'f'csf ) 

TECX-tOO-empdt 

c 

return 

end 
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Appendix  F:  Output  Parameters  for  PLOT.TMP 


1 

u 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

1  :  External  reactivity  {absolute  value) 

2  :  Effective  reactivity  insertion  ( absolute  value  ) 

3  :  Feedback  reactivity  (absolute  value) 

4  :  Reactor  thermal  power 

5  :  Coolant  temperature  at  core  inlet  (A') 

6  :  Coolant  temperature  at  core  outlet  (A') 

7  :  Coolant  average  temperature  in  the  core  (A') 

S  :  Fuel  average  temperature  in  the  core  (A) 

0  :  Cladding  average  temperature  (A'l 

10  :  Pressure  loss  in  the  core  annulus  (:Y /  m~) 

11  :  Pressure  loss  in  the  core  (A :(m~) 

12  :  Total  core  pressure  loss  (JV  jm2) 

13  :  TE  hot  shoe  temperature  in  the  TEM  pump  (A') 

14  :  TE  cold  shoe  temperature  in  the  TEM  pump  (A') 

15  :  Electric  current  in  the  TEM  pump  (amp) 

16  :  Power  consumed  in  the  TEM  pump  (&W') 

17  :  Total  pressure  loss  in  the  pipes  (.V /m~) 

18  :  Pressure  loss  in  the  heat  exchanger  ( N/m 

19  :  Total  pressure  losses  in  the  system  (N/m2) 

20  :  Pump  head  losses  in  the  system  (N/m2) 

21  :  Primary  mass  flow  rate  (kg /sec) 

22  :  Input  power  to  the  main  TE  generators  (k\V) 

23  :  Electric  power  output  (fcW) 

24  :  Rejection  power  through  the  main  radiator  (k\V) 

25  :  Conversion  Efficiency  (°c) 

26  :  Electric  current  for  a  single  TE  generator  f  amp ) 

27  :  Electric  voltage  for  a  single  TE  generator  (volt ) 

28  :  Hot  shoe  temperature  of  the  main  TE  converter  (A') 

29  :  Cold  shoe  temperature  of  the  main  TE  converter  (A) 

30  :  TEM  Pump  efficiency  (cr) 

31  :  Secondary  coolant  inlet  temperature  (A') 

32  :  Secondary  coolant  outlet  temperature  (A) 

33  :  Secondary  mass  flow  rate  (kg  / sec ) 

34  :  Average  temperature  of  the  main  radiator  (A') 
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Appendix  G:  Transient  Simulation  for  Negative  Reactivity  Insertion 
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Figure  G.l.  Transient  Case  with  Negative  Reactivity  Insertion 
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Figure  G.3.  Temperature  at  Key  Points  as  a  Function 
Transient  Case  with  Negative  Reactmu 
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